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Abstract 




The  effects  of  iii5)urities  - hydrogen,  water  vapour  and  sodium  chloride  - 
on  shock  wave  stability  and  structure  in  ionizing- argon  and  krypton  flows  were 
investigated  by  using  a dual-wavelength  Mach-Zehnder  interferometer  in  conjunction 
with  the  OTIAS  10  cm  x l8  cm  (4'"^x  7"T  Hypervelocity  Shock  Tube. 

At  shock  Mach  number  of  about  15,  when  pure  argon  or  krypton  is  used 
a.s  a test  gas,  the  resxilting  translational  shock  front  develops  sinusoidal 
oscillations.  The  subsequent  ionization  relaxation  region,  the  electron- cascade 
front  and  the  quasi-eqx;illibrium  state  also  exhibit  unstable  characteristics.  The 
addition  of  small  amounts  of  hydrogen  (^0.5^  by  pressiore)  to  the  test  gas 
stabilizes  the  entire  flow  and  at  the  same  time  drastically  shortens  the  relaxa- 
tion length  to  about  l/3  of  the  original  value  in  argon  and  2/3  in  krypton.  The 
addition  of  about  1^  of  water  vapour  (by  pressure)  to  the  argon  test  gas  also 
stabilizes  the  flow  and  shortens  the  relaxation  region.  A thin  coating  of 
dissolved  sodium  chloride  in  water  on  the  shock  ttibe  wall,  which  on  evacuation 
leaves  a rough  surface  of  crystals,  however,  does  not  seem  to  have  any  effect  on 
the  flow. 


Although  the  small  in5)urities  of  hydrogen  and  water  vapour  have  a 
significantVffect  in  stabilizing  the  shock,  the  overall  values  of  density  and 
electron  number  density  in  the  quasi- equilibrium  region  remain  nearly  xhe  same. 


The  present  experimental  results  provided  the  data  for  a theoretical 
study  which  took  Into  account  atoni-atom  collisions,  electron-atom  collisions  and 
radiation-energy  losses  to  deduce  a precise  value  for  the  atom-atom  collision 
excitation  cross-section  constants  for  axgon  and  krypton  having  values  of 
^Ar-Ar  “ ^ 10“ ^9  cn^/eV  and  = 1.2  x 10“19  cm^/eV,  respectively. 


"Kr-Kr 
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Notation 


a 


At 


c 

e 


g 

g 

h 

k 

Kr 

Kr"* 

L 

m 

a 

m 

e 

M 


n 


Si 

^a 


Atom 

Argon  atom 
Argon  ion 
Speed  of  light 
Electron 

Atomic  statistical  weight 
Gaunt  factor 
Planck  constant 
Boltzmann  constant 
Krypton  atom 
Krypton  ion 

Distance  travelled  by  light  source  in  shock  tube 

Mass  of  atom 

Mass  of  electron 

Shock  wave  Mach  nimiber 

Refractive  index 

Nimiber  density  of  atoms 

Number  density  of  electrons 

Electron  number  density  production  rate 

Pressure 

Rate  of  radiant  energy  loss/unit  volume  by  line  radiation 

Rate  of  radiant  energy  loss/unit  volvaae  by  continuum  radiation 

Rate  of  total  radiant  energy  loss 

Electron  atom  momentum  transfer  cross-section 

Electron  ion  momentum  transfer  cross-section 

Nondimensionalized  fringe  shift 
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Notation  - Continued 


a- a 


e-a 


T 

T. 


e 

u 

Z 


eff 


a 


0 

0 


P 

u 


(0 


p 


Atom-atom  excitation  collision  cross  section  constant 
Electron-atom  inelastic  collision  cross  section  constant 
Ten?)erature 

Heavy  particle  temperatiire 
Electron  temperature 
Particle  velocity 
Relaxation  length 
Effective  nuclear  charge 
Degree  of  ionization 

Energy  transfer  rate  due  to  inelastic  collision 
Energy  transfer  rate  due  to  elastic  collision 
Ionization  teii5>erature 
Cutroff  frequency 
Density 

Angular  frequency 
Plasma  frequency 


Siib  scripts 
a Atom 

e Electron 

i Ion 

s Shock,  plasma  species 
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1.  IHTRQDUCTION 


The  ionization  of  monatomic  gases  behind  strong  shock  waves  is 
now  understood  quite  well  owing  to  considerable  theoretical  and  experimental 
work  devoted  to  the  subject.  Aft»:r  some  detailed  theoretical  and  experimental 
considerations,  Petschek  and  Byron  (Ref.  1)  concluded  that  the  mechanisms 
leading  to  equilibrium  ionization  were  first  controlled  by  inpurities  and 
then  dominated  by  electron-atom  collisions.  Harwell  and  Jahn's  (Ref.  2) 
experimental  work  supported  the  proposal  of  Bond  (Ref.  3)  that  the  initial 
ionization  proceeded  via  a two-step  atom-atom  collision  process;  that  is, 
electronic  excitation  to  the  first  excited  state  followed  by  ionization. 
Subsequent  investigators  (Refs.  4-8)  all  adopted  the  model  for  achieving 
ionization  relaxation  as  first  due  to  atom-atom  and  then  electron-atom 
collisions.  However,  there  are  qiaite  a nimiber  of  discrepancies  in  the  value 
used  for  the  argon  atom-atom  excitational  collision  cross-section  constant 
(Refs.  2,  4,  9,  10). 

Harwell  and  Jahn  (Ref.  2)  observed  that  vmknown  inpurities  in 
the  test  gases  had  profound  effects  on  the  ionization  rates  at  concentration 
as  low  as  a few  parts  per  million.  Morgan  and  Morrison  (Ref.  4)  using  a 
relatively  sinple  theory  and  based  on  some  characteristic  excitation  energy 
tried  to  fit  Harwell  and  Jahn's  (Ref.  2)  experimental  resiiLt.  They  hypo- 
thesised that  sodium  and  nitrogen  inpurities  had  a definite  effect  whereas 
hydrogen  and  oxygen  had  no  significant  effect  on  the  rate  of  ionization. 

These  are  somewhat  contrary  to  the  findings  of  Bristow  (Ref.  ll),  who  showed 
that  oxygen  had  no  observable  effect  but  hydrogen  did  have  a significant 
effect  in  reducing  the  ionization  relaxation  time.  Schneider  and  Park  (Ref. 
12)  studied  the  ionization  rates  using  NaCl- contaminated  argon  as  a test  gas. 
Their  measurements  showed  that  an  increase  in  the  rate  of  ionization  was  due 
to  electron  detachment  of  negative  chlorine  ions  produced  from  the  sodiian 
chloride  vapour.  They  also  concluded  that  atmospheric  inpiurities  such  as 
H2O,  N2  and  CO2  had  insignificant  effect  on  the  rate  of  ionization.  This 
result  is  contradicted  by  the  present  experiments. 

Enomoto  (Ref.  8)  included  the  pertxirbations  caused  by  the 
laminar  boundary  layer  in  his  work  and  showed  that  the  ionization  relaxation 
time  shortened  as  the  shock  tube  diameter  decreased.  (He  also  showed  that 
the  ionization  relaxation  time  was  dependent  on  the  initial  pressiore  for  a 
given- Mach  nunber  and  shock  tilie  diameter.)  McLaren  and  Hobson  (Ref.  lO) 
also  included  boundeiry  layer  effects  in  their  theoretical  considerations. 

It  should  be  noted  that  the  above  authors  all  adopted  Mirels'  (Ref.  13) 
theory  which  considered  only  ideal  gases. 

D'yakov  (Ref.  l4)  initiated  the  study  of  shock  wave  stability 
and  later  ^wan  et  al  (Ref.  15)  and  Fowles  (Ref.  16)  re-examined  the  subject 
and  came  up  with  a condition  for  stability  for  real  gases  at  very  high  Mach 
number  when  the  Hugoniot  changes  slope  and  curvature.  Griffiths  et  al  (Ref. 
17)  have  also  looked  into  this  problem.  Levine  (Ref.  18)  studied  shock  wave 
stability  due  to  turbulent  mixing  at  the  contact  surface.  The  above  investi- 
gators were  not  involved  with  the  same  type  of  instability  first  observed  by 
Bristow. 

The  present  research  was  designed  to; 

(1)  Experimentally  Investigate  the  effects  of  inpurities  on  shock  wave 
stability. 
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(2)  Provide  experimental  data  for  the  detennination  of  atom-atom 
collision  excitation  cross-section  constants,  and 

(3)  Provide  experimental  data  for  the  choosing  of  a rediation-loss 
model. 


2.  SHOCK  STRUCTURE  ANALYSIS 


Before  considering  the  shock  structure  equations,  the  colli- 
sional  ionization  processes  employed  in  this  analysis  will  be  discussed. 
Immediately  behind  the  Rankine-Hugoniot  shock  front,  the  gas  is  heated  to 
high  temperature  by  a few  translational  collisions,  which  increase  the 
kinetic  energy  of  the  atoms.  Then  via  a two  step  atom-atom  collisional 
process,  which  takes  some  atoms  to  their  first  excited  state,  and  siibsequent 
collisions  ionize  the  atoms.  In  this  manner  the  first  few  electrons  are 
produced.  The  next  phase  is  dominated  by  the  much  more  efficient  electron- 
atom  collisional  process,  leading  to  a very  rapid  production  of  electrons. 

As  the  electron  number  density  reaches  a certain  level,  the  three  body 
(electron-ion-electron)  recombination  process  becomes  important  and  this 
is  followed  by  a quasi- eqviilibrium  state.  Subsequently,  the  electron 
number  density  drops  due  to  radiant  energy  losses  as  electrons  and  atoms 
recosmbine  and  de-excite.  In  the  radiative  cooling  region,  the  plasma  is 
almost  in  an  equilibrium  state.  It  can  be  concluded  that  the  entire  shock- 
wave  structure  can  be  divided  into  two  regions,  the  ionization  relaxation 
region  and  the  radiative  cooling  region.  The  former  region  is  dominated 
by  collisional  processes  while  the  latter  region  is  dominated  by  radlatlbh 
energy  losses. 


For  a two- component  mixtiire,  i.e.,  argon  with  hydrogen  or 
krypton  with  hydrogen,  the  following  collisional  ionization  processes  are 
considered: 

a + a -a  + a + e 

a + e a + e + e 

a + H a^  + H + e 

H + a - + a + e 

H + H + H + e 

H + e + e + e 


Here,  we  assume  that  the  ionization  of  hydrogen  molecviles  starts  effectively 
behind  the  strong  shock  since  the  relaxation  length  for  dissociation  of 
hydrogen  molecules  is  very  short  compared  with  that  of  ionization  (Ref.  19). 
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The  following  assumptions  are  made  in  the  theoretical  analysis; 

(l)  the  flow  is  steady  and  quasi-one-dimensional  when  viewed  in  shock-fixed 
coordinates;  (2)  the  gas  is  singly  ionized;  (3)  the  heavy  particles  (atoms 
and  ions)  have  the  same  temperatures  and  velocities;  (4)  the  heavy  particles 
and  electrons  have  Maxwellian  velocity  distributions  corresponding  to  the 
local  atom  temperature  and  electron  tenperature;  (5)  the  electron  diffusive 
velocity  with  respect  to  the  heavy  particle  velocity  is  very  small  and  can 
be  neglected;  (6)  the  effects  of  electric  and  magnetic  fields  on  the  shock- 
wave  structure  can  be  neglected;  and  (7)  the  effects  of  viscosity  and  thermal 
conductivity  axe  very  small  and  can  be  neglected. 

Under  these  assumptions,  the  governing  equations  which  are  applicable 
to  the  nonequilibrivun  relaxing  flow  behind  the  shock  wave  are  (Ref.  20): 


= 0 

(2.1) 

P)  = 0 

(2.2) 

= -Qip 

(2.3) 

II 

a 

Ve 

(2.4) 

dT 

IT  ^ ^ ^ ^ 

(2.5) 

where  Qr  is  the  rate  of  radiation  energy  loss  of  the  plasma,  n©  is  the 
electron  ntmiber  density  production  rate,  T|  is  the  energy  transfer  rate  given 
to  electrons  by  elastic  collision,  0 is  the  energy  transfer  rate  given  to 
electrons  by  inelastic  collisions,  p is  the  plasma  density,  u is  the  plasma 
velocity,  P is  the  flow  pressure,  H is  the  total  enthalpy,  a is  the  degree 
of  ionization  [a  = ne/(ne  + na) ] , ng  and  Te  are  the  electron  nimriber  density 
and  tenperature,  respectively. 


Expauiding  the  above  equations  for  a two- component  mixture  of  plasma, 
they  may  be  written  as  (Ref.  21) 


( V ^ = ^01  V 


^“^a2  ^ ^e2^'^  ~ ^02  ^^02 


) -I  m .(n  . + n .)u^  + n .kT^  + (n  . + n . 
1 aj'  aj  ej'  ej  e ' aj  ej 

3=1  2 

‘ I { ”aj  "OJ  “OJ  * "OJ  } 


)kT, 


(2.6) 

(2.7) 


(2.8) 
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life 

j=l 


1 “^5 

7?  m .(n  . + n .)vr  + ^ (n  . + n .)Kr  u + 

2 aj'  aj  ev  2 ' aj  ej''  a 


tL 

^n.KE  u + n.ke.u  1=-)  (Or,.) 

2 ej  e ej  j jj 


J=1 


to  = A 


CL  / \ • 

to  (”e2^)  = 


o dl  j r~> 

2”e’*Sr*'‘e“eS'^  <"3  * '>3> 

J=1 

where  subscripts  1 and  2 denote  components  of  the  gas  mixture. 


el 

‘e2 

2 


(2.9) 

(2.10) 
(2.11) 

(2.12) 


Since  the  percentage  of  the  hydrggen  inpurity  is  small,  the  kinetic 
energy  and  momentum  transfers  between  the  two  types  of  heavy  particles  are 
neglected. 


The  rate  of  radiant  energy  loss,  Qr,  consists  of  the  rates  of  energy 
loss  by  continuum  radiation,  Qc,  and  by  line  radiation,  Qb*  continuum 

radiation  in  a plasma  is  caused  by  recombination  processes  (free-bound  radia- 
tion) and  free-free  transitions  and  is  given  by  (Ref.  22) 


e e 


(2-13) 


where  Zeff  is  the  effective  nuclear  charge.  Meiners  and  Weiss  (Ref.  23)  found 
that  Z^ff  = 1.67  for  argon,  and  z|ff  = I.65  for  krypton. 


The  energy  loss  by  line  radiation  due  to  bound-bound  transitions 
can  be  written  as  (Ref.  24) 


(2.14) 

k 

where  is  the  energy  level  of  the  k'*'^  excited  state  and  k_i  the  radiative 
transition  rate  coefficients  and  na(k)  the  number  density  of  atoms  at  the  k^^ 
excited  state  and  is  given  by 


^ Z'  h^  ^3/2  2 (ke-Ej^)/kT^ 

' ' -n  ^ 


2gj  27r  m^KT^  ; 


n e 
e 


(2.15) 


where  gj^  and  g];  are  the  atomic  statistical  weights  of  the  k^^  excitation 
level  and  of  the  ion,  respectively. 
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For  krypton,  we  assumed  that  the  rate  of  radiant  energy  loss  by 
line  radiation  is  equal  to  that  for  continuum  radiation.  The  same  assun5)tion 
is  also  made  for  the  hydrogen  plasma. 


The  elastic  energy  transfer  rate,  t),  is  the  siim  of  the  rates  of 
thermal  energy  given  to  the  free  electrons  by  electron-atom  collisions,  i]ea> 
and  by  electron-ion  collisions,  T)ei- 


Tiea  ^ei  related  to  their  corresponding  momentum  transfer 
cross-sections  Qea  and.  Qei  (Ref.  21) 


■ea 


= 6 n n 


e a V tt 


~ k(T  - T ) 
m 'a  e' 
a 


Ti  . = 6 n 
/ei  e 


m KE 
e e 


TT  J 


^ k(T  - T ) 
m 'a  e' 
a 


(2.16) 

(2.17) 


The  inelastic  energy  transfer  rate,  0,  is  the  sum  of  the  rates  of 
thermal  energy  given  to  the  free  electrons  by  electron- atom  and  electron- 
ion-electron  inelastic  collisions  0ea  and  by  bremsstralilung  (^. 

0g^  and  0p  are  given  by  (Refs.  21,22) 


0 = - (n  ) 

ea  ' e'  e 


ke  + ^ kT^ 


(2.18) 


Stt  n ^ e^  kT  / Qir  m n1/2 

^ m 2 h V ^e  / 
e 

Further  discussions  on  the  theoretical  considerations 
21  and  25 . 


(2.19) 

may  be  found  in  Refs. 


3-  EXPERIMENTAL  FACILITIES  AMD  TECHNIQUES 

Shock  ttibes  have  been  used  very  successfully  in  the  study  of  plasma 
flows  (Refs.  1,  2,  5}  10,  11).  Diagnostic  techniques  are  quite  different 
and  can  be  divided  into  optical  interferometry  (Refs.  5>  11) > electrostatic 
probes  (Refs.  1,  10),  microwaves  (Refs.  2,  9)  and  photoelectric  interferometry 
(Ref.  26). 


The  present  experiments  were  conducted  in  the  UTIAS  combustion 
driven  10  cm  x l8  cm  Hypervelocity  Shock  Tube.  A Mach-Zehnder  interferometer 
with  a 23-cm  diameter  field  of  view,  using  a pulsed  laser  light  soiurce  was 
used  as  the  main  diagnostic  equipment.  Since  the  original  construction  of 
the  facility  (Refs.  27,  28),  several  modifications  were  made  which  substan- 
tially iii5)roved  the  performance  and  ease  of  handling.  The  details  are 
described  in  Refs.  11  and  29. 
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3.1  Shock  Tube 


The  shock  tube  can  be  divided  mainly  into  four  sections:  driver 
section,  channel  or  driven  section,  test  section  and  dump  tank.  The  driver 
section  has  a length  of  1.42m  and  an  inner  diameter  of  15.24  cm.  A O.38  mm 
diameter  iiqpulsively  heated  tungsten  wire  is  located  at  the  centre  of  the 
driver.  An  experimentally  determined  stoichiometric  mixture  of  7.5^  (by 
pressure)  oxygen,  1%  hydrogen  diluted  by  72. 5?^  heliian  and  ^ hydrogen  was 
used  as  a driver  gas.  When  a 13  kV  potential  was  applied  impulsively  across 
the  tungsten  wire,  the  gas  mixture  was  ignited  by  combustion.  The  final 
driver  pressure,  P4,  was  increased  to  about  6.6  times  that  of  the  initial 
gas  mixture.  A typical  pressure  trace  of  the  driver  is  shown  in  Fig.  1. 

It  can  be  seen  clearly  how  the  pressure  is  built  15)  after  combustion  and 
drops  suddenly  eifter  the  stainless  steel  diaphragm  natures. 

The  ratio  of  driver  pressure,  P4,  to  initial  argon  test  gas 
pressure.  Pi,  was  found  experimentally  by  Bristow  (Ref.  U)  to  be  related 
to  the  incident  shock  Mach  nxmiber  by  the  empirical  expression. 


in 


= (0.422M  + 4.02),  for  10  < M <24,  argon 

S “ S “* 


For  the  krypton  test  gas,  there  were  not  a sufficient  number  of 
runs  to  formiiLate  a simiilar  expression  over  a range  of  Mach  nimnbers.  But 
it  is  found  that  for  Mach  nurdber  Ms  ~ 15)  the  following  compression  holds. 


in 


(0.422M  + 2.92),  for  M ~ 15)  krypton 
s s 


The  diaphragms  were  made  of  annealed  stainless  steel  type  SS304 
with  an  ultimate  tensile  strength  of  85)000  psi.  They  were  scribed  according 
to  the  following  equation  (Ref.  U), 


where  cr  is  the  ultimiate  tensile  strength, 
t is  the  diaphragm  thickness, 
d is  the  unsupported  diameter,  = 18.I  cm, 
h is  the  thickness  after  scribing, 

K and  n are  values  determined  experimentally,  to  ensure  proper 
petallng  after  btirst. 

A transition  section  is  located  between  the  driver  and  the  driven 
section  where  the  cross-section  is  changed  from  the  circular  15.24  cm 
diameter  to  the  rectangular  10  cm  x I8  cm  cross-section. 

The  driven  section  has  a length  of  J3.2m  with  ports  located  1.22m 
apart  which  can  be  used  to  house  gauges  to  monitor  the  passage  of  the  shock 
wave  or  to  record  its  pressure  history.  The  10  cm  x I8  cm  test  section  is 
1.22m  long  eind  Incorporates  two  Interferometer- quality  glass  windows  with 
a 20  cm  diameter  aperture  and  9.5  cm  thick.  They  are  polished  to  a ssx>othness 
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of  l/4  wavelength  at  6328^  with  wedge  aJigles  not  exceeding  1 sec  of  arc. 
Cross-hairs  8,2  cm  apart  axe  movmted  on  the  windows.  A dunp  tank  located 
at  the  end  of  the  channel  serves  as  a safety  feature  in  preventing  a 
reflected  shock  wave . 

The  vacuum-punp  system  of  the  shock  tube  consists  of  a Kinney 
type  KS-47  mechanical  pun?),  a Hereaus  Rootes-type  VP-RG-350A  pun?)  and  a 
CVC  type  PMCS-6B  diffusion  plan?),  which  are  capable  of  bringing  the  pressure 
down  to  3 X 10"5  torr  with  an  outgassing  rate  of  about  2 x 10"5  torr/min. 
Prolonged  pun?)down  period  does  not  seem  to  in?)rove  the  ultimate  vacuum. 

A mass  spectrometric  study  (Ref.  29)  of  the  present  vacuum  system  showed 
that  the  main  desorption  constituents  were  water  vapour,  air  and  some 
hydrocarbons.  There  were  no  detectable  amounts  of  sodium  chloride,  as 
suggested  by  Schneider  and  Park  (Ref.  13). 

3.2  Interferometer 


The  interferometer  used  is  of  the  Mach-Zehnder  type  (Ref.  28) 
and  has  a 23- cm  diameter  field  of  view.  A TRG  104A  rviby  laser,  pxilsed  by 
a Pockels  cell  Q-switch  and  fitted  with  a second  harmonic  generator, 
provided  the  694 3A  and  3471 •5A  con?)onent  wavelength  light  source.  A 
schematic  of  the  interferometer  system  is  shown  in  Fig.  2. 

A detailed  description  of  the  construction,  modification  and 
operation  of  the  interferometer  can  be  found  in  Refs.  28,  11,  29.  The 
basic  working  principle  of  the  interferometer  is  outlined  below. 

Light  shines  through  the  input  mirror  (Fig.  2)  which  is  so  placed 
as  to  reflect  the  light  towards  the  collimating  parabolic  mirror  1.  It  is 
then  divided  into  two  portions  at  splitter  1,  one  travels  throvigh  the  com- 
pensation chamber  and  onto  mirror  2,  while  the  other  via  mirror  1 passes 
through  the  test  section.  The  two  beams  then  rejoin  at  splitter  2 and 
form  interference  patterns  depending  on  the  optical  path  lengths  each  hsis 
travelled.  Parabolic  mirror  2 then  sends  the  image  through  the  output 
mirror  to  the  camera.  Inside  the  camera  are  located  a 60-40  beam  splitter, 
a plane  mirror  and  two  line-filters  to  separate  out  the  caii?)onent  wavelengths. 
The  film  used  was  Kodak  Royal  X Pan  (1250  ASA) . DK-50  proved  to  be  the 

proper  developer  with  a developing  time  of  about  8 minutes. 

The  attachment  of  the  ruby  laser  directly  onto  the  frame  of  the 
interferometer  greatly  reduced  the  difficulty  of  adjustment  (Ref.  29).  As 
has  been  noticed  before  (Ref.  11),  air  currents  and  teinperature  variations 
have  significant  effects  on  the  orientation  of  the  fringes.  Although  the 
interferometer  is  covered  by  canvas  and  foam  sponge,  it  is  still  subject 
to  ten?)erature  gradients.  The  teinperature  rise  is  attributed  to  the  heat 
generated  by  the  surrounding  electronic  equipment  and  the  white  light  source 
used  in  setting  and  focussing  the  central  order  fringes.  In  order  to  min- 
imize the  movements  of  the  fringes,  the  use  of  the  white  light  source  is 
restricted  and  the  air- conditioner  is  shut  off  to  avoid  air  currents  prior 
to  each  run.  The  ultimate  way  to  ensure  stable  fringe  settings  might  be 
an  insulated  housing  for  the  interferometer. 

3 *3  Velocity  Measurement 

Accurate  measurement  of  the  shock  velocity  is  necessary  in  order 
to  provide  a precise  shock  Mach  nuntoer  for  the  theoretical  calculations. 
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As  vd.ll  be  shovm  later,  Mach  number  has  a significant  effect  in  determining 
the  atom-atom  collision  excitation  cross-section  constant.  A schematic 
circuit  diagram  for  velocity  measurement  is  illustrated  in  Fig.  3. 

As  mentioned  before  (Sec.  3-l)j  there  are  ports  located  along  the 
channel  and  in  the  test  section  of  the  tube.  These  ports  axe  used  to  hovise 
piezoelectric  pressure  transducers  (Atlantic  Research  type  LD-25)  to  monitor 
the  arrival  of  the  incident  shock  wave.  Their  outputs  are  aii?)lified  and 
used  to  start  and  stop  electronic  counters  thereby  giving  several  velocity 
measurements  of  the  shock  as  it  travels  along  the  channel.  From  the  velocity 
measurements  it  is  possible  to  determine  the  shock-wave  attenuation.  At  a 
shock  Mach  number  Mg  ~ l6  travelling  into  argon  at  5 torr,  the  attenuation 
is  about  0.7^  per  meter  while  at  Mg  ~ 15  in  5 torr  krypton,  it  is  about  1^ 
per  meter.  That  is,  the  attenuation  of  the  shock  wave  is  more  severe  at 
the  lower  Mach  number  runs  in  krypton  than  at  the  higher  Mach  number  runs 
in  argon.  A detailed  boundaxy  layer  analysis  would  probably  yield  the 
reasons  to  explain  these  facts. 

Even  though  the  shock  wave  velocity  history  along  the  channel 
was  knovm,  the  velocity  measurement  taken  across  stations  I and  J (Fig.  3) 
was  used  to  determine  the  incident  shock  Mach  nunber  at  the  test  section. 

This  was  done  for  consistency.  A Kistler  pressvjxe  transducer  vms  used  at 
station  M,  right  at  the  centre  of  the  test  section  vd.ndows,  to  record  the 
pressure  history  of  the  shock  induced  flow.  Typical  pressure  traces  for 
argon  and  krypton  runs  are  shovm  in  Fig.  4.  Figure  4a  shows  the  pressure 
history  of  a shock  wave  Mg  = 15*9  travelling  into  5 *14  torr  of  argon.  The 
Inward  shift  of  the  trace  indicates  the  arrival  of  the  shock  wave.  It  is 
not  possible  to  pinpoint  the  exact  location  of  the  ionization  shock  front 
from  the  pressure  trace  as  the  ionization  relaxation  time  is  only  about 
4 fisec  which  is  well  within  the  rise  time  of  about  10  nsec  for  the  pressure 
transducer.  The  symbol  S represents  the  positions  of  both  the  translation 
and  ionization  shock  fronts.  From  the  trace  the  pressure  behind  the  shock 
front  has  an  average  value  of  about  1700  torr  over  the  flow  period.  The 
predicted  value  by  assuming  the  flow  in  ionization  equilibrium  vrithout 
considering  radiation  losses  is  l88o  torr.  The  total  test  time  available 
is  estimated  to  be  about  200  nsec  for  this  case  including  real  gas  and 
viscoios  effects  but  not  the  turbvilent  diffusion  of  the  contact  front. 
Experiments  show  that  the  actual  test  time  for  the  conditions  studied  is 
only  about  100  nsec.  It  is  very  difficult  to  accurately  locate  the  position 
of  the  contact  front  on  the  pressure  trace  as  the  pressure  across  the  contact 
front  is  constant.  Nevertheless  the  beginnings  of  the  contact  front  at  100 
nsec  axe  so  indicated.  Figure  4b  shows  the  pressure  history  of  a shock 
wave  at  Mg  = 15.0  travelling  into  5.15  torr  of  krypton.  A pressure  jurp  of 
about  280  times  is  expected  and  shows  up  accordingly  in  the  trace.  The 
fluctuations  of  the  filtered  pressure  trace  probably  come  about  mainly 
from  the  natural  frequencies  associated  with  the  transducer  diaphragm  and 
cannot  be  explained  in  detail  as  relating  to  the  flow. 

3-4  Admission  of  Test  Gases  and  lapurltles 

The  test  gases  used  vrere  of  the  highest  purity  level  available 
commercially.  The  argon  test  gas  had  a quoted  purity  level  of  99.9998^t> 
vdiereas  that  of  the  krypton  test  gas  was  99 *995^  • The  higher  grade  argon 
did  not  have  any  noticeable  effects  in  reducing  shock  wave  instability. 
Neither  did  the  admission  of  the  test  gas  through  ports  located  upstream 
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or  downstream  of  the  tube  or  purging  the  whole  tube  with  test  gas  before 
runs  have  any  influence  on  stability.  The  initial  pressure  of  the  test  gas 
was  monitored  by  a Wallace  & Tieman  gauge  (0-50  torr  range)  and  a more 
accurate  measurement  of  the  initial  pressure  was  made  using  an  oil  manometer. 

For  the  runs  that  involved  the  addition  of  hydrogen  as  an  impurity, 
the  technique  used  by  Brimelow  (Ref.  29)  was  employed.  The  same  technique, 
however,  proved  unsuitable  for  admitting  water  vapoxir. 

It  was  found  that  the  water  when  exposed  to  low  pressure  will 
evaporate,  cool  and  freeze.  The  ice  thus  formed  will  stay  in  the  solid 
state  for  a long  time  under  low  pressure.  A simple  device  was  used  later 
to  introduce  water  vapour.  A measured  amount  of  distilled  water  was  put 
inside  a saiall  metal  capsule  which  was  then  lowered  into  liquid  nitrogen. 

The  vrater  was  frozen  almost  instantly.  The  capsule  was  then  connected  to 
an  inlet  with  a valve  near  the  evacuated  test  section  (station  J,  Fig.  3). 

The  valve  on  the  inlet  was  cpened  and  since  the  volume  of  the  capsule  was 
very  small  compared  with  the  tube,  the  pressure  inside  was  down  to  that 
of  the  tube  in  a very  short  time.  The  water  inside  the  capsiile  was  still 
frozen  and  presimiably  very  little  of  it  was  lost  by  evaporation.  The  valve 
was  then  shut  and  the  capsule  heated  to  vaporize  the  water.  At  the  same 
time  the  pumps  were  shut  off  and  part  of  the  test  gas  was  let  into  the  tube 
and  mo^i  tored  by  a McLeod  gaiJge.  The  valve  was  opened  again  to  let  in 
the  water  vapour.  Some  pressure  measxirements  from  the  McLeod  gauge  are 
given  in  Table  1. 

There  are  a few  things  worth  noting  from  these  measurements. 

First,  not  all  the  water  inside  the  capsiale  was  converted  to  water  vapour. 

For  example,  if  0.5  ml  of  water  were  all  vaporized,  then  we  would  have  a 
measurement  of  SOOp  instead  of  250p.  Second,  as  can  be  seen  from  Table  1, 
all  the  readings  are  decreasing  with  time.  This  probably  indicates  that 
part  of  the  water  vapour  introduced  was  being  adsorbed  by  the  shock  tube 
walls  forming  mo3  ecilLar  layers  which  <jo)jld  not  be  jneasured  by  the  pressure 
gauge.  The  reason  for  introducing  part  of  the  test  gas  first  and  then  water 
vapour  was  to  saturate  the  wall  surfaces  with  the  test  gas  thereby  hopefully 
allowing  all  the  water  vapour  to  be  part  of  the  test  gas.  As  can  be  seen, 
this  was  not  the  case.  Apparently,  water  vapour  has  a tendency  of  displacing 
all  other  gas  molecules  that  cling  to  the  wall  sxurface.  The  previously 
mentioned  miass  spectromietric  study  seems  to  s\pport  the  above  concept. 

Since  there  is  no  accurate  means  in  the  present  instrumentation 
equipment  of  measuring  the  exact  amount  of  water  vapour  in  the  test  gas 
mixtiore,  the  last  readings  in  Table  1 were  used.  This  is  by  no  means 
precise  for  when  the  balance  of  the  test  gas  was  introduced,  a purging 
effect  coxild  take  place  and  possibly  some  of  the  adsorbed  water  vapour 
evolved.  However,  as  the  present  experiments  were  meant  to  be  qualitative, 
the  readings  chosen  will  probably  sezire  the  purpose. 

For  one  experiment,  a thin  coating  of  sodium  chloride  crystals 
was  grown  on  the  bottom  side  of  the  tiibe  in  the  test  section  before  the  run. 

A normal  saline  solution  was  applied  inside  the  tube  and  then  the  pressure 
inside  was  lowered  slowly  by  controlled  vacuum  pumping.  As  the  water 
vaporized,  a thin  and  uniform  coating  was  formed  on  the  test- section  lower 
wall  (Fig.  5).  The  results  will  be  discussed  in  the  subsequent  sections. 
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4.  MALYSIS  OF  IKTEEFEROGRAMS 


Basically,  an  interferogram  consists  of  Bright  and  dark  fringes 
as  shown  in  Fig.  6.  Sequential  numbers  are  assigned  to  the  fringes  as 
shown.  It  gives  a relative  measurement  of  refractive  index  between  a 
reference  and  an  unknown  region.  The  (x,y)-plane  consists  of  a pre-shock 
region,  a translational  shock  front  S,  an  ionization  relaxation  region  Xe, 
an  electron  cascade  front  E and  a post- shock  region.  The  overall  shock  wave 
SE  is  moving  to  the  left  into  the  gas  at  rest. 


31) 


The  refractive  index  of  a gas  is  related  to  its  density  by  (Ref. 


(n  - 1) 

s 


(4.1) 


vdiere  n is  the  refractive  index, 

is  the  Gladstone-Dale  constant  for  species  ’s', 

p**  is  the  corresponding  density  for  species  's', 
s 


and  subscript  s denotes  the  coniJonents  of  the  gas.  in  the  case  of  a singly 
ionized  monatomic  gas  that  consists  of  neutral  atoms,  ions  and  electrons. 
The  contributions  of  the  electronically  excited  atoms  to  the  refractive 
index  are  not  included  (Ref.  11).  In  the  rims  that  involve  hydrogen  as  a 
small  iii5)urity,  its  contributions  to  the  refractive  index  were  neglected. 


As  can  be  seen  from  the  pre- shock  region  in  Fig.  6,  the  density 
there  is  constant  so  there  is  no  change  in  refractive  index  and  the  fringes 
show  i^)  as  uniformly  horizontal.  The  almost  discontinuous  downward  shift 
of  fringes  signifies  the  location  of  the  incident  translational  shock  front. 
This  is  followed  by  a region  where  atom-atcjm  collisions  take  place.  As 
this  process  is  not  very  efficient  in  producing  electrons  (i.e.,  the  flow 
properties  are  change  very  slowly),  the  change  of  refractive  index  in 
this  region  is  very  slow  and  as  a result  the  fringes  gradually  shift 
upweird.  As  more  and  more  electrons  are  formed,  they  become  very  effective 
collision  partners  and  a rapid  cascading  of  electrons  occurs.  This  results 
in  a rapid  change  in  the  refractive  index  and  the  fringes  correspondingly 
rise  suddenly. 


Fringe  shifts  make  it  possible  to  measure  the  total  plasma 
density  p and  the  electron  (ion)  number  density  ne.  The  fringe  shifts 
are  defined  as  the  difference  of  fringe  number  (or  orders  of  interference) 
at  a particular  location  where  the  original  fringe  is  being  displaced  by 
a new  fringe.  For  exaaple  if  the  location  of  point  X (Fig.  6)  was  originally 
occupied  by  fring  No.  15  and  now  that  location  registers  fringe  No.  13,  then 
at  that  location,  there  is  a fringe  shift  of  2.  The  nondimensional  (normal- 
ized by  the  fringe  spacing  5 to  6 or  6 to  7,  say)  fringe  shift,  S,  is  related 
to  the  change  of  refractive  index  by  (Ref.  29) 


® ■ i;  ("j  ■ “i' 

(4.2) 

S "I  [(Hj  - 1)  - (n^  - D] 

(4.3) 
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where  L is  the  distance  (IO.I6  cm)  the  light  source  travelled  across  the 
test  section, 

X is  the  wavelength  of  the  lig^it  source  expressed  in  cm  and  subscripts 
i,  J denote  the  reference  and  investigated  regions,  respectively. 

Expanded  for  the  case  of  a singly  ionized  plasma,  Eq.  4.1  becomes 
(Ref.  32)  2 

(n  - 1)  = p[K^(l  - a)  + (4.4) 

2u) 


vdiere  - Gladstone-Dale  constant  for  atoms, 

Kj  - Gladstone-Dale  constant  for  ions, 
a - degree  of  ionization, 

Wp  - plasma  frequency, 
and  u)  - angular  frequency. 

Equation  4.4  can  also  be  expressed  as 


(n 


1)  = pK^  - pa 


where  m is  the 
a 

After 


mass  of  atom. 

some  algebraic  manipulation 


of  Eqs.  4.3  and  4.5, 


Pj  = "i 


. ^1^2^  - "2"l\ 
LK 


and 


n . = n . + 
ej  ei 


V2^  - Ka2^i\ 


LKm 


where  n - electron  nunfcer  density. 


“1  ■ (“a  - Vl  * 


X 10"^^ 


m 


°2  ■ - %>2  ♦ 


U.>(85  X 10"^** 


m 


= "l^A2  - °2Ka1 


(4.5) 


(4.6) 

(4.7) 


and  subscripts  loSnd  2 correspond  to  wavelengths  Xq  (6943.0  x 10"®  cm)  and 
>2  (3471.5  X 10"®  cm)  respectively.  The  values  of  the  constants  for  argon 
and  krypton  in  these  equations  are  listed  in  Table  2. 
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Ideally,  if  we  know  the  density  and  electron  number  density  in 
the  reference  region,  and  the  Gladstone-Dale  constants  of  all  the  constituents 
in  the  unknown  region,  then  by  measuring  fringe  shifts  at  points  of  interest, 
we  can  calculate  the  density  and  electron  number  density  at  those  points. 
However,  due  to  optical  imperfections  there  may  be  some  fictitious  fringe 
shifts.  The  no-flow  interf erograms , which  are  taken  immadiately  before 
firing,  together  with  the  flow  interf erograms,  will  provide  the  true  fringe 
shifts,  and  thus  the  actual  changes  in  flow  properties  can  be  calculated. 

4.1  Gladstone-Pale  Constants 

As  noted  in  Sec.  4,  a knowledge  of  the  Gladstone-Dale  constants 
for  all  the  species  in  the  plasma  is  essential  in  determining  the  changes 
in  flow  properties.  Peck  and  Fisher  (Ref.  33)  used  a Michelson  interfero- 
meter to  measure  the  refractivity  of  argon  and  came  up  with  a dispersion 
formula  for  argon.  Chashchina  and  Shreider  (Ref.  34)  used  the  method  of 
spectral  line  shift  to  determine  index  of  refraction  for  krypton  and  a 
dispersion  formula  was  also  derived.  Using  Cauchy's  dispersion  equation 
(Ref.  35)  for  krypton,  the  refractive  index  obtained  was  slightly  different, 
but  the  overall  results  in  density  and  electron  number  density  evaluated  are 
almost  the  same.  Expressions  for  the  free-electron  refractive  index  can  be 
found  in  Refs.  11  and  36.  Following  the  expression  for  free  electrons  in 
Ref.  36,  a constant  of  4.485  x 10“14  ^2  obtained  and  used  in  Sec.  4, 

Eq.  4.5.  This  value  differs  from  that  given  in  Ref.  36  but  agrees  with 
Ref.  11. 


The  values  for  the  argon  and  krypton  ions,  however,  are  not  that 
well  agreed  upon.  Alpher  and  White  (Ref.  36)  used  a Slater  screening- 
constant  theory  and  calculated  the  refractive  index  of  the  argon  ion  to  be 
0.67  times  that  of  the  neutral  atom..  Later  investigators  (Refs.  11,26,37) 
all  obtained  a value  of  0.72  using  the  same  theory.  Bristow  (Ref.  ll)  did 
an  experimental  determination  of  the  polarizability  for  singly  ionized 
argon  and  obtained  a ratio  of  0 . 65  for  the  argon  ion  to  the  argon  atom. 
However,  the  value  obtained  by  Alpher  and  White  is  still  being  widely  used 
(Refs.  26,30,38).  Igra  (Ref.  38)  compared  the  resulting  density  profiles 
by  using  various  ratios  from  0.25  to  1.0  and  concluded  that  the  difference 
is  very  small. 

Since  there  is  no  experimental  result  available  for  the  krypton 
ion  and  using  Igra's  corparison,  the  ratio  of  the  index  of  refraction  of 
the  krypton  ion  to  its  atom  was  taken  to  be  0.75,  as  calculated  by  using 
Slater  screening- constant  theory.  It  should  be  noted  that  Meiners  and 
Weiss  (Ref.  26)  used  the  ratio  of  1 in  their  calculation. 

4.2  Digitizing  Equipment 

The  present  scheme  of  analyzing  interferograms  was  developed  by 
Whitten  (Ref.  29).  It  requires  a transformation  of  the  fringe  locations 
into  digital  form  with  respect  to  a chosen  reference  point.  For  exaaple, 
if  we  choose  a particular  point  in  an  interferogram  as  the  reference  point 
and  assign  it  to  be  the  origin  of  an  (x,y) -coordinate  system,  then  the 
fringes  can  be  thought  of  as  lines  or  loci  of  points.  Figxire  7a  shows  an 
interferogram  of  a shock  wave  travelling  from  right  to  left.  The  p re- shock 
region,  the  translational  shock  front,  the  ionization  relaxation  region. 
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the  electron  casceuie  front  and  the  radiative  cooling  region  can  al l be 
seen  clearly.  Cross-hairs  are  located’ 8.2  cm  apart  and  can  be  used  as 
reference  points.  The  cross-hair  at  the  lower  right  hand  corner  of  this 
interferogram  was  chosen  as  one  of  the  reference  points  and  also  as  the 
origin  of  an  (x,y)- coordinate  system.  The  corresponding  digitized  version 
of  the  pre-shock  region,  the  translational  shock  front,  the  ionization 
relaxation  region,  the  electron  cascade  front,  the  radiative  cooling  region 
and  the  reference  point  (origin)  are  shown  in  Fig.  7b.  With  all  four  inter- 
ferograms  (2  no- flow  and  2 flow)  of  one  experiment  digitized  in  this  fashion 
with  respect  to  the  same  reference  point  (origin),  we  can  determine  the 
actual  fringe  shift  when  the  same  (x,y) -coordinate  of  any  point  in  the  flow 
and  no- flow  pictiires  is  cco^iared. 

To  transform  fringes  into  digital  form,  a digitizer  was  used. 
First,  it  involves  picking  a common  reference  point  in  all  four  interfero- 
grams  and  in  this  case,  the  cross-hairs  were  used.  Then  a hand-held  ctirsor 
is  lased  to  trace  the  fringes,  and  their  locations  with  respect  to  the 
reference  point  are  recorded  directly  on  a magnetic  tape.  The  recorded 
data  are  later  transferred  from  tape  onto  cards.  An  IBM  l800  Plotter  is 
then  used  to  reproduce  the  (igitized  version  of  the  interferograms . There 
are  of  co\arse  some  errors  in  tracing  the  exact  locations  of  fringes.  By 
con{)aring  the  digitized  and  original  interferograms , we  can  pick  out  the 
points  that  are  not  aligned  with  the  original.  Since  the  points  are  in 
digital  form  on  cards  already,  they  can  be  corrected  and  the  error 
reduced  until  a good  reproduction  is  obtained.  The  original  interfero- 
grams are  enlarged  to  about  twice  their  actual  size  which  makes  the  whole 
digitizing  and  correction  procedure  much  easier. 

The  whole  set  of  points  are  then  linearly  interpolated  into 
uniform  steps.  Fringe  nujribers  (orders)  are  also  interpolated  and  assigned 
to  each  (x,y)-coordinate . With  the  flow  and  no- flow  pictures  at  the  two 
wavelengths  now  digitized  and  fringe  nxmibers  assigned,  fringe  shifts  at 
any  point  can  be  calculated  and  in  turn  the  changes  in  the  flow  properties 
at  that  point.  The  entire  flow  region  is  then  known. 


5 . RESULTS  AND  DISCUSSIONS 


Bristow  (Ref.  ll)  first  observed  that  in  high  Mach  nuniber  (Mg  ^ 
15)  argon  rims,  a ccm5)lex  system  of  flow  distturbances  persist.  Figure  8 
shows  two  interferograms  taken  simultaneously  using  two  wavelengths, 

6943. OA  and  3471. 5A.  An  incident  shock,  Mach  number  Mg  = l4.7,  is  moving 
from  right  to  left  into  argon  at  an  initial  pressure  Pi  = 4.08  torr  and 
teii5)erature  Ti  = 298. 0 K.  Clearly  seen  axe  the  sinusoidal-type  oscillations 
of  the  translational  shock  front  S,  where  the  fringes  change  abruptly, 
followed  by  a relaxation  region,  terminating  in  the  electron- cascade  front 
E,  followed  by  disturbances  in  the  fringes  at  this  front  and  beyond  into 
the  q\iasi-equilibrium  region.  Subsequent  investigators  (Refs.  29,30>38) 
using  the  same  experimental  facilities  all  reported  similar  findings. 

Figure  9 shows  an  interferogram  of  a shock  wave  at  Mg  = 15*2  travelling 
from  right  to  left  into  the  krypton  test  gas  initially  at  Pi  = 5-07  torr 
and  Ti  = 295-4  K.  Evidently,  the  same  type  of  disturbances  also  exist  here. 
Other  investigators  (Refs.  22,36)  who  used  a rotating  mirror  camera  to 
record  time-resolved  dxial-wavelength  interferograms  of  the  ionized  argon 
flow,  did  not  come  across  this  phenomenon.  The  major  apparent  difference 
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was  in  the  cross-sections  of  the  shock  tubes.  Alpher  and  White  (Ref.  36)  used 
a tiibe  8.25  cm  x 8.25  cm,  while  Oettinger  and  Bershader  (Ref.  22)  used  a tube 
5 cm  X 5 cm.  The  UTIAS  tube  has  a cross-section  of  10  cm  x I8  cm. 

Morse  and  Ingrad  (Ref.  39)  have  treated  the  problem  of  ar^difi cation 
of  acoustic  disturbances  in  a colli si on- dominated  plasma  where  initially  the 
electron  temperatxire  Te  is  much  higher  than  the  heavy  particle  temperature  Ta- 
in our  case,  immediately  behind  the  translational  shock  front  in  ionizing 
krypton  and  argon,  the  reverse  is  true  (see  Fig.  10)  but  the  ten^ierature 
difference  between  electrons  and  heavy  particles  is  large.  Based  on  this  fact, 

Bristow  postulated  that  such  a condition  may  be  favourable  to  generating  ion- 
acoustic  disturbances.  Since  the  calculations  by  Belozerov  and  Measures  (Ref. 

4o)  for  strong  shocks  into  hydrogen  have  shown  that  the  electron  temperature  is 
essentially  the  same  as  the  heavy  particle  temperature  in  the  post- shock  relaxa- 
tion region,  Bristow  then  added  a small  percentage  of  hydrogen  to  the  argon  test 
gas  for  the  purpose  of  reducing  the  temperature  difference  between  the  electrons 
and  the  heavy  particles.  As  will  be  shown  later,  the  addition  of  hydrogen  had 
a very  significant  effect. 

But  the  results  obtained  by  Whitten  (Ref.  29)  and  Brimelow  (Ref.  30) 
at  lower  Mach  numiber  argon  rtins  (Ms  " 13),  disagreed  with  Bristow's  postulate. 

Figure  11  shows  two  simultaneous  interferograms  of  a shock  wave  (Ms  = 13-0) 
travelling  from  right  to  left  abound  a sharp  corner  in  ionizing  argon  with  an 
initial  pressure  and  temperature  of  5-01  torr  and  296.8  K,  respectively.  Dis- 
regarding the  effects  of  the  expansion  wave  at  the  corner,  we  can  see  that  in 
the  free  stream,  the  translational  shock  front  is  planar,  the  siibsequent 
ionization  relaxation  region  and  the  electron  cascade  front  are  all  stable. 

Figure  12  shows  the  result  at  a higher  Mach  numiber  (Mg  = 16.I)  argon  run.  A 
comparison  of  Figs.  11,  8 and  12  reveals  that  as  the  incident  shock  Mach  numiber 
increases,  or  the  degree  of  ionization  increases,  the  more  severe  become  the 
disturbances.  This  strongly  suggests  that  the  disturbances  observed  are 
associated  with  significant  ionization.  Figure  13c  shows  a strong  shock  wave 
(Ms  = 16.6)  travelling  into  nitrogen  at  an  initial  pressure  Pi  = 2.26  torr. 

Nitrogen  is  dissociated  but  not  ionized  in  this  case.  As  can  be  seen,  the  flow 
is  stable  and  free  of  any  distvirbances . As  further  proof  that  shock  wave 
disturbances  sere  associated  with  significant  ionization  only,  series  of  runs 
were  carried  out  using  krypton  as  the  test  gas.  Krypton  can  be  readily  ionized 
and,  as  can  be  seen  from  Fig.  9,  the  dlstiorbances  are  produced  as  expected.  The 
above  mentioned  experiments  were  all  conducted  in  the  same  manner  using  the 
same  facility. 

In  Ref.  21,  there  are  some  detailed  discussions  on  the  factors 
affecting  the  theoretically  calculated  relaxation  length,  Xg  (that  is,  the 
distance  between  the  shock  front  S to  the  point  of  maximvim  electron  number 
density  E,  see  Fig.  6).  It  is  shown  in  Ref.  21  that  for  a fixed  shock  Mach 
numiber  Mg  and  an  initial  pressure  Pi,  Xg  theory  = Xg  (S^.aj  Sg.g,  Q^a,,  Qeij 
Qg) . That  is,  the  relaxation  length  depends  on  the  atom-atom  excitation 
collision  cross-section  constant,  the  electron  atom  inelastic  collision  cross- 
section  constant,  the  electron  atom  momentm  transfer  cross-section,  the  electron 
ion  miomientum  transfer  cross-section  and  the  lOfjradiatipn  energy  lossi  A.  study  of 
their  relative  importance  shows  that  is  the  main  and  most  sensitive  factor 
in  calculating  Xg. 


l4 


Based  on  the  experimental  data  from  pure  test-gas  runs  and  with 
hydrogen  inmurities,  best-fit  values  were  obtained  of  Sa-a  for  argon, 

= 1.0  X 10“^9  cm2/eV,  and  for  krypton,  = 1.2  x 10"l9  cm^/eV.  The 

available  information  in  the  relatively  long-duration  qiiasi-equilibrium  region 
helped  in  the  choice  of  the  radiation  model.  The  initial  conditions  for  the 
cases  studied  are  given  in  Table  3*  Some  representative  results  are  given 
below.  Additional  results  can  be  found  in  Appendix  A. 


5.1  Pure  Argon 

Figure  l4  shows  the  interferogram  for  Case  1,  Table  3}  and  Fig.  15 
shows  the  corresponding  plot  of  the  nondimensional  electron  number  density 
ne/ne  E through  the  shock  wave  transition  and  beyond,  as  well  as  the  nondimen- 
sional density  p/pE»  As  can  be  seen  in  Fig.  l4,  the  translational  shock  front 
is  sinusoidal  and  disturbances  show  up  as  wavy  fringes.  The  lower  half  of  the 
shock  front  appears  imusually  thick  owing  to  nonplanarity.  It  substantiates 
the  hypothesis  that  the  shock  wave  is  not  one-dimensional  and  may  even  be  three- 
dimensional  in  nature.  As  the  translational  shock  fronts  always  show  up  sinu- 
soidal in  shape  (Figs.  8,  9,  12,  13  and  l4),  it  is  very  possible  that  the  fronts 
are  also  wavy  in  a direction  normal  to  the  optical  viewing  path.  To  avoid  non- 
planar  effects  on  the  shock  front,  measurements  were  taken  in  the  most  quiescent 
regions  which  usually  also  avoided  the  major  oscillations  throi;igh  the  relaxation 
region  and  beyond. 

The  post  shock  values  of  ne/ng  ^ in  Fig.  15  are  also  oscillatory  and 
reflect  the  waviness  of  the  fringes  in  Fig.  l4,  as  expected.  The  interfero- 
metric data  points  were  taken  1 mm  apart  for  convenience  and  more  or  fewer  points 
can  be  evaluated  depending  on  the  work  required. 

The  position  X = 0 (Fig.  15)  is  where  the  translational  shock  front  is 
located.  The  position  X = Xe,  where  ng/ng  g = 1 is  the  position  of  the  electron 
cascade  front  or  the  end  of  the  transition  region.  Here  the  flow  quantities 
achieve  their  maximum  equilibrium  values.  Radiation  losses  are  insignificant 
in  this  region  but  beyond  this  point  radiation  losses  became  iuportant  and  the 
electron  number  density  falls.  The  fits  to  the  experimental  data  of  ng  and  p 
axe  quite  reasonable  despite  their  oscillations . The  experimental  relaxation 
length  is  2.0  cm  conpared  with  the  theoretically  determined  value  of  2.1  cm  when 
the  value  of  is  chosen  to  be  1.0  x 10“^9  cD^/eV.  Figure  l6  shows  a plot 

of  the  nondimensional  degree  of  ionization  cc/ct£  with  distance  X.  The  post  shock 
values  after  3 cm  are  in  paxticilLarly  good  agreement  with  the  analysis.  The 
effects  of  radiation  losses  or  radiative  cooling  are  very  much  in  evidence  and 
lend  sipport  to  the  choice  of  the  radiation-loss  model. 


The  fact  that  an  error  of  ± 0.3  in  the  measured  Mach  number  was  possible 
(as  discussed  in  Sec.  5*5)  leads  to  the  result  that  the  determined  value  of  S^r-Ar 
will  change  slightly.  The  measured  Mach  number  for  this  case  is  Mg  = 15.9  (solid 
line).  Two  additional  cxirves  are  plotted  for  Mg  = l6.2  (dashed)  and  15.6  (dot- 
dashed)  . In  order  to  obtain  a best  fit  to  the  interferometric  data,  the  corre- 
sponding values  have  to  be  changed  to  0.9  x 10“^9  cm^/eV  and  1.1  x 10"19 

cn^/eV,  respectively.  Since  the  equilibrium  (real  gas  without  radiation)  values 
for  Mg  =15.9  were  used  for  normalizing,  the  three  cxzrves  appear  separated.  Only 
this  case  was  used  to  demonstrate  the  effects  of  a possible  error  in  Mach  number 
on  the  effective  value  of  S^_^.  The  error  bars  shown  in  Figs.  15,  I6  and  in 
subsequent  cases  are  due  to  rringe  location  reading  errors  only  as  will  be  discussed 


15 


later  in  Sec.  5.5*  Table  4 gives  the  eqitLlibrium  values  of  pg, 
for  the  cases  considered. 

Case  2 (Figs.  17,  l8)  has  similar  initial  conditions  to  Case  1 (see 
Table  3)  and, as  expected,  the  resiilts  are  also  similar.  The  experimentally 
measured  relaxation  length  (1.9  cm)  conpares  favourably  with  the  theoretically 
determined  value  (2.0  cm).  The  oscillations  of  data  points,  especially  in  the 
post- shock  region,  are  very  much  in  evidence  due  to  wavy  fringes  in  the  inter- 
ferograms.  The  measured  maxlmijm  electron  number  density  (degree  of  ionisation) 
overshoots  the  theoretical  value,  but  in  general,  the  agreement  is  reasonable. 
Case  3 (Figs.  19,  20)  has  a slightly  higher  Mach  number  and  as  a result  the 
relaxation  length  is  shorter  (1.7  cm  for  theory  and  1.8  cm  for  experiment). 
Despite  the  oscillations,  the  agreement  between  theory  and  experiment  is  quite 
good.  It  should  be  noted  that  although  the  agreement  between  the  experimental 
results  and  the  analysis  for  the  above  cases  is  fair  to  good,  it  is  not  so  for 
Case  4 (Figs.  21,  22)  and  Case  5 (Figs.  23,  24).  Figures  21  and  22  show  the 
resvilts  for  a shock  Mach  nimiber  Ms  = 13.0.  It  can  be  seen  that  the  theoretical 
and  interferometrically  measured  relaxation  length  differ  by  2.0  cm  over  a length 
of  about  11  cm.  One  probable  reason  for  this  is  that  at  lower  Mach  nimiber  the 
changes  of  density  and  electron  number  density  are  relatively  small  and  as  a 
result  the  changes  in  refractive  index  become  the  same  order  as  the  errors 
associated  with  the  measurements.  Figures  23  and  24  (Case  5)  show  the  results 
of  a higher  Mach  number  run  (Mg  = 17-2) . The  agreement  in  the  quasi- equilibrium 
region  is  quite  reasonable,  but  not  for  the  relaxation  length.  The  maximum 
electron  number  density  measured  and  calculated  are  also  in  poor  agreement. 
Enomoto  (Ref.  8)  pointed  out  that  for  higher  Mach  number  runs,  the  relaxation 
time  (length)  is  significantly  shortened  when  wall -boundary-layer  effects  are 
included.  The  lower  initial  pressure  may  also  contribute  to  the  shortening  of 
relaxation  length.  As  mentioned  before,  the  main  impurities  in  the  desorbed 
gases  consist  of  water  vapour,  air  and  some  hydrocarbons.  The  lower  the  initial 
pressure,  the  higher  the  percentage  of  water  vapoiir  in  the  test  gas.  As  will 
be  discussed  later,  water  vapour  has  a significant  effect  in  reducing  the  relaxa- 
tion length.  This  probably  helps  to  explain  the  discrepancy  in  the  calculated 
and  measured  relaxation  length.  However,  the  reason  for  the  lajrge  discrepancy 
in  the  maximvm  electron  number  density  for  this  case  is  still  unknovm.  It  may 
be  necessary  to  fit  a new  value  of  the  atom-atom  collision  excitation  cross- 
section  constant  for  higher  (and  lower)  shock  Mach  number.  Consequently,  until 
the  values  of  the  quantities  used  to  determine  Xe  (see  Sec.  5)  are  accurately 
known,  we  cannot  claim  that  is  a universal  constant. 

Table  5 lists  the  predicted  and  measured  relaxation  lengths  for  the 
cases  studied.  The  conparisons  between  theory  and  experiment  are  quite  good 
for  Cases  1,  2 and  3-  For  Case  4 where  the  shock  Mach  number  is  much  lower 
than  the  others,  the  agreement  is  not  so  good.  Case  5 has  a lower  initial 
pressure  and  a higher  Mach  number.  The  difference  between  predicted  and  measured 
relaxation  length  is  quite  severe.  The  effects  of  wall-boundary- layer  suggested 
by  Enomoto  (Ref.  8)  and  the  possibility  of  higher  inpurity  level  in  the  test  gas 
probably  cause  the  discrepancy. 

5 .2  P’jre  Krypton 

Experimental  resvilts  for  krypton  are  quite  sparse.  Hopefully,  the 
present  data  can  provide  additional  valuable  information.  Figure  25  shows  an 
interferogram  of  a shock  wave  moving  into  5.07  torr  of  krypton  at  a Mach  nunber 
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^ =15-2  (Case  6).  The  sinusoidal  shock  front  is  very  much  in  evidence  and 
so  is  the  varying  shock-front  thickness  which  further  substantiates  that  the 
wave  is  in  fact  nonplanar.  The  distiirbances  in  the  electron  cascade  front  are 
strikingly  similar  to  those  in  pure  argon  runs. 

The  closeness  of  the  electron  cascade  front  to  the  translational  shock 
wave  at  the  wall  is  also  worth  noting.  Figures  26  and  27  show  the  plots  of  the 
nonnaliaed  electron  number  density,  density  and  degree  of  ionisation  profiles. 
The  agreement  between  the  interferometrically  measured  results  and  analysis  is 
very  good  for  the  chosen  value  of  the  atom-atom  collisional  excitation  cross- 
section  constant  = 1.2  x 10-19  cn^/eV.  The  oscillation  of  the  data 

points  correspond  with  the  wavy  fringes  due  to  disturbances  and  not  due  to 
errors  in  the  data  reduction.  Figures  28  and  29  show  the  results  for  another 
krypton  case.  The  length  of  the  quasi- equilibrium  region  provides  information 
on  the  choice  of  the  radiation-loss  model  used  in  the  theoretical  calculation. 

A conparison  of  the  measured  and  calculated  relaxation  lengths  can  be  found  in 
Table  5«  As  can  be  seen,  the  agreement  between  theory  and  experiment  for  both 
cases  are  very  good. 

5 .3  Inpurities 

Inpiorities  have  long  been  inportant  factors  in  the  study  of  shock- 
wave  structure  and  stability.  Petschek  and  Byron  (Ref.  1)  and  Harwell  and  Jahn 
(Ref.  2)  all  observed  that  impurities  have  a profound  effect  on  the  rate  of 
ionization  and  probably  provide  the  initial  mechanism  leading  to  equilibrium 
ionization.  However,  little  agreement  exists  on  which  impurities  have  or  have 
not  an  effect  on  the  ionization  rate. 

Morgan  and  Morrison  (Ref.  k)  using  a relatively  simple  theory  and 
based  on  some  characteristic  excitation  energy  tried  to  fit  Harwell  and  Jahn's 
(Ref.  2)  experimental  results.  They  hypothesized  that  a nitrogen  impurity  had 
a definite  effect  whereas  oxygen  and  hydrogen  impurities  had  no  significant 
effects  on  the  rate  of  ionization.  The  results  shown  here  (from  Ref.  ll) 
contradict  part  of  their  hypothesis.  Figure  13a  shows  an  interferogram  of  a 
shock  wave  travelling  into  argon  at  an  initial  pressure  of  2.85  torr  with  an 
oxygen  impurity  (Q.46^  by  pressxire) . The  disturbances  are  very  much  in 
evidence  and  the  interferogram  resembles  a pure-gas  rxm  (Figs.  8,  9,  12  6ind  l4). 
A comparison  of  the  mieasured  and  calculated  (assuming  pirre  argon  as  the  test 
gas)  relaxation  lengths  shows  the  agreement  to  be  reasonable.  This  proves  that 
an  oxygen  impurity  has  no  significant  effect  on  shock- wave  structure  and  stabi- 
lity and  agrees  with  Morgan  and  Morrison's  hypothesis.  Figure  13b  shows  an 
interferogram  of  a Mach  17.2  shock  wave  travelling  into  argon  with  a helium 
impurity.  The  results  are  also  simdlar  to  those  for  pure  argon  runs  and  that  a 
helium  impurity  has  no  significant  effect  on  either  shock  wave  structure  or 
stability.  Figure  13c  shows  a strong  shock  wave  travelling  into  nitrogen.  The 
test  gas  was  dissociated  but  not  ionized.  It  can  be  seen  clearly  that  the  flow 
behind  the  shock  front  is  very  stable  as  expected.  No  experiment  was  done  using 
nitrogen  as  an  impurity  in  argon  as  a test  gas,  owing  to  the  negative  results 
from  Fig.  13g.  It  is  very  doubtfiil  that  nitrogen  would  have  any  significant 
effect  on  argon  shock- wave  structure  and  stability.  An  additional  test  would 
clarify  this  point. 

Schneider  and  Park  (Ref.  12)  studied  experimentally  the  ionization 
rates  using  argon  as  a test  gas  with  sodium  chloride  as  an  impurity.  Their 
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measurements  showed  that  the  increase  in  ionization  rate  was  due  to  electron 
detachment  of  negative  chlorine  ions  produced  from  the  sodium  chloride  vapoizr. 

The  present  stucJy  shows  that  a hydrogen  impurity,  which  can  be  readily  ionized 
to  produce  electrons,  increases  the  ionization  rate  of  argon.  Both  studies 
attribute  the  increase  in  the  ionisation  rate  to  the  electrons.  It  has  beer- 
shown  that  electrons  are  very  efficient  collision  partners  which  can  lead  to  a 
rapid  production  of  new  electrons  until  equilibrium  is  achieved.  Although  it 
is  well  known  that  sodium  chloride  from  the  atmosphere  can  cling  to  the  shock- 
tube  walls  during  a routine  cleanup,  oixr  mass  spectrometric  study  failed  to 
reveal  any  appreciable  amounts  of  soditim  chloride  in  the  desorbed  gas.  As 
sodium  chloride  can  be  easily  ionised  to  produce  electrons,  one  experiment  was 
performed  by  having  sodiimn  chloride  crystals  on  the  shock- tube  wall.  A normal 
saline  solution  was  brushed  onto  the  shock- tiibe  wall  and  then  by  carefully 
monitoring  the  vacuum  pumps,  the  pressiore  inside  the  tube  was  lowered  slowly. 

The  water  evaporated  and  a rough  coating  of  crystals  was  formed  on  the  wall 
surface.  It  was  hoped  that  when  the  shock  wave  travelled  over  the  crystals 
it  woxjld  scrub,  heat  and  ionise  the  salt  which  might  affect  the  shock  structure. 
However,  upon  examination  of  the  resulting  interferogram  (Fig.  30) j it  did  not 
exhibit  features  different  from  a normal  piu:’e  argon  run. 

5-3-1  Hydrogen  Impurity 

Bristow  (Ref.  ll)  first  discussed  the  use  of  hydrogen  addition  with 
a view  to  eliminate  the  inexplicable  flow  distiirbances  both  within  and  behind 
the  shock-wave  structure.  Figiore  31  shows  five  interferograms  of  similar  shock 
strength  and  initial  conditions.  They  only  differ  in  the  amount  of  hydrogen 
addition  to  the  test  gas.  Figure  31a  shows  the  run  with  an  almost  pure  argon 
test  gas;  the  rest  all  have  different  amounts  of  hydrogen  addition.  As  the 
amount  of  hydrogen  added  to  the  test  gas  increases,  the  flow  disturbances 
diminish  and  the  relaxation  length  shrinks.  It  was  fo\mid  by  Bristow  (Ref.  ll) 
that  0.4^  (by  pressiire)  of  hydrogen  addition  to  the  argon  test  gas  can  completely 
eliminate  the  disturbances.  Figures  31d  and  31e,  which  represent  two  runs  with 
0.%  and  lio  hydrogen  addition,  respectively,  show  a straight  translational  shock 
front  along  with  a stable  relaxation  region,  electron  cascade  front  and  radiative 
cooling  region  compared  to  Fig.  31aj  which  has  almiost  no  hydrogen  addition. 

Besides  stabilizing  the  flow,  hydrogen  also  shortens  the  entire  shock  transition. 
Figure  31  also  illustrates  how  the  relaxation  length  decreases  as  the  amount  of 
hydrogen  addition  increases.  It  was  found  that  by  adding  a small  amouint  of 
hydrogen  (0.4^  by  pressiore)  as  an  impurity  to  the  argon  test  gas,  the  relaxation 
length  was  reduced  to  l/3  of  the  pure-gas  value.  This  result  is  quite  different 
from  the  argon-xenon  mixtures  used  by  Kelly  (Ref.  9).  He  showed  that  the  addi- 
tion of  0.1^  or  0.48^  of  xenon  to  the  argon  test  gas  did  not  substantially 
change  the  relaxation  length  from  that  for  pure  argon.  The  reason  is  that  the 
collision  excitation  cross-section  constant  for  xenon  is  much  smaller  than  for 
argon,  while  the  constant  for  hydrogen  is  much  larger.  In  addition,  the  mass 
of  the  hydrogen  atom  is  markedly  smialler  (l/4o)  than  that  of  argon  atom  and  the 
ionization  energy  for  atomic  hydrogen  is  lower  than  that  for  argon.  Tnese  three 
characteristics  of  a hydrogen  impurity  in  an  argon  test  gas  probably  account  for 
the  large  reduction  in  the  relaxation  length. 

Figure  32  shows  a similar  comparison  of  the  effects  of  hydrogen  addition 
to  krypton  as  a test  gas.  It  can  be  seen  that  with  0.4^  hydrogen  impurity  in 
the  krypton  test  gas  (Fig.  32b),  the  flow  is  the  most  stable  one.  Although 
hydrogen  has  a similar  effect  to  argon  in  reducing  the  relaxation  length  of 
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ionizing  krypton,  it  does  not  completely  eliminate  the  flow  disturbances. 
Furthermore,  with  more  than  0.k%  hydrogen  added,  it  seems  to  enhance  or  even 
produce  other  types  of  disturbances  (Figs.  32c  and  32d) . 

The  reason  why  an  addition  of  hydrogen  can  stabilize  the  flow  in 
argon  is  still  unknown  at  this  time.  Although  pure  krypton  rims  exhibit 
similar  traits  as  for  pure  argon  runs,  the  effects  of  hydrogen  addition  are  not 
identical  as  new  types  of  instabilities  appear.  The  reason  why  hydrogen  behaves 
more  favourably  with  argon  but  not  with  krypton  is  not  known.  It  should  be 
pointed  out  that  although  a .smaJ  1 addition  of  a hydrogen  inpurity  in  an  argon 
or  krypton  test  gas  significantly  reduces  the  relaxation  length  and  eliminates 
the  distxirbances,  it  does  not  affect  the  overall  flow  properties  in  the  quasi- 
equLlibrivon  region. 

Calculations  were  made  for  the  cases  of  Ms  = l6.2,  Pp  = 5 >17  torr  and 
Ms  =17.6,  Pp  = 3.12  torr  (Cases  8 and  9)  and  with  a sman  amount  of  hydrogen 
(0.4^  by  pressure)  as  an  impurity  in  argon.  The  resulting  analytical  and  experi- 
mental density,  electron  nuiriber  density  and  degree  of  ionization  variations 
through  the  ionization  relaxation  and  radiative  cooling  regions  are  shown  in 
Figs.  33-36.  The  shock  transition  values  are  particularly  good  for  the  stronger 
shock  case  (Figs.  35  and  36).  On  the  whole,  the  agreement  between  the  measured 
and  predicted  relaxation  lengths  (see  Table  5)  and  the  quasi-equilibrium  post- 
shock value  are  qxute  satisfactory.  Better  agreement  mi^t  have  been  obtained 
by  including  the  appropriate  values  of  the  refractive  indices  for  the  hydrogen 
impurity  in  the  interferometric  equations  (Eqs.  4.6  and  4.7).  The  assimption  of 
equal  velocities  and  temperatures  for  the  heavy  particles  of  argon  and  hydrogen 
in  the  theory  may  also  contribute  some  error. 

Similar  calc\ilations  were  also  made  for  the  cases  (Cases  10  and  ll) 
of  krypton  with  0.38^  (by  pressure)  of  added  hydrogen  impurity.  The  compari- 
sons of  measured  and  predicted  density  and  electron  number  density  from  the 
pre-shock  region  through  the  ionization  relaxation  and  radiation  cooling  regions 
are  shown  in  Fig.  37-  The  ionization  relaxation  length  for  this  case  is  8 mm 
when  measured  experimentally  and  the  theory  predicted  9 nnn  (Table  5)  • Both  the 
measured  density  and  electron  number  density  profiles  are  below  the  predicted 
values.  Since  the  degree  of  ionization  profile  is  a function  of  both  density 
and  electron  numiber  density  profile,  the  measured  values  overshoot  the  predicted 
ones  (Fig.  38).  The  agreement  between  experiment  and  theory  is  much  better  in 
Case  11.  The  maeasured  density  profile  in  Fig.  39  agrees  very  well  with  the 
predicted  values.  The  agreement  in  the  long  quasi-equilibrium  region  is 
especially  good.  The  interferometrically  determined  electron  numiber  density 
profile  falls  sli^tly  short  of  the  predicted  values,  but  the  overall  trend  is  in 
reasonable  agreemient.  Figure  4o  shows  the  corresponding  degree  of  ionization 
profile  and  as  can  be  seen,  the  agreement  between  theory  and  experiment  is  also 
very  good,  confinnj.ng  the  choice  of  the  theoretical  model. 

5.3.2  Water  Vapour  Impurity 

The  possibility  of  impurities,  especially  water  vapour,  being  present 
on  any  shock  tube  surface  cannot  be  disregarded.  Even  thou^  the  wall  was 
cleaned  as  thorovighly  as  possible  before  each  run,  the  presence  of  water  vapour 
moleciiles,  typically  at  a density  of  10^  cm”^  at  a pressure  of  10“°  torr,  are 
always  there  as  shown  by  our  mass  spectroscopic  study.  Schneider  and  Park  (Ref. 
13)  pointed  out  that  water  vapour  had  no  significant  effect  on  the  rate  of  ioni- 
zation of  argon.  There  is,  however,  no  mention  about  what  other  effects  water 
vapour  may  have. 
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The  use  of  water  vapour  as  an  impurity  was  originally  intended  to 
study  wall  effects  in  the  ionization  relaxation  region  (See.  The  method 

used  for  admitting  hydrogen  inpurity  was  first  tried.  It  involved  using  a 
precalibrated,  known  Wallace-Tieman  gaiage  volume  reservoir  to  determine  the 
amount  of  water  (by  pressure)  being  admitted  into  the  test  gas.  It  was  later 
found  that  most  of  the  water  vapour  remained  inside  the  reservoir.  Even  the 
vapour  in  the  test  gas  probably  formed  molecular  layers  of  water  on  the 
surfaces,  which  could  not  be  detected  by  a pressure  gauge.  Consequently, 
it  was  not  possible  to  know  the  exact  amount  of  water  vapour  added  to  the 
test  gas . Figure  4l  shows  a resulting  interferogram  when  an  unknown  but 
"small"  amount  was  presented  in  the  argon  test  gas.  The  unstable  shock  front 
is  moving  to  the  left,  leaving  behind  a non-uniform  region.  Such  disturbances 
were  net  observed  before.  The  same  experiment  was  repeated  using  the  reservoir 
but  no  water  was  added.  The  argon  test  gas  was  admitted  into  the  tube  via  the 
reservoir.  This  procedure  had  not  been  xased  before  as  the  test  gas  was  usually 
admitted  through  another  inlet.  The  resulting  interferogram  is  shown  on  Fig. 
k2.  A conparison  of  Figs.  4l  and  k2  shows  an  almost  identical  disturbance 
even  though  the  latter  one  was  intended  to  be  a pure  argon  run.  Part  of  the 
no- flow  region  is  missing  and  the  reason  for  its  blockage  is  not  clear.  It 
was  concluded  that  perhaps  the  adsorbed  or  remaining  water  inside  the  gauge 
reservoir  from  the  previous  run  must  have  been  driven  into  the  tube  when 
admitting  the  argon  test  gas  for  the  second  run.  When  the  original  method  of 
admitting  the  argon  through  another  inlet  was  used  again,  then  the  usual 
interferogram  (Fig.  12)  was  again  obtained.  Additional  experiments  using  water 
vapoxar  as  an  irpurity  are  still  being  carried  out  at  the  present  time  (Ref.  ^l) 
Preliminary  results  confim  that  very  small  amounts  of  water  vapour  added  to 
the  argon  test  gas  will  produce  disturbances  similar  to  those  shown  in  Figs. 

4l  and  42. 


A more  efficient  method  of  admitting  water  vapour  was  later  devel- 
oped (Sec.  3*4).  Although  it  is  still  not  possible  to  make  an  accurate  measure 
ment  of  water  vapoxir  concentration  in  the  test  gas,  a qualitative  measurement 
can  be  obtained.  The  difficulty  is  mainly  due  to  the  fact  that  water  vapour 
tends  to  adsorb  on  the  surfaces  and  no  siople  measuring  device  is  available 
that  can  measure  the  low  concentration  of  water  vapour  (~  0.1'^  of  5 torr  by 
pressure)  of  interest  in  the  present  work. 

Figure  43  shows  an  interferogram  of  a perfectly  stable  shock  wave 
in  argon.  A small  amount  (0.25  ml)  of  distilled  water  was  introduced  into 
the  shock  tube.  If  all  the  water  turned  into  vapour  then  the  additional 
expected  pressure  would  be  300p.  However,  30  seconds  after  the  water  vapour 
was  admitted  it  indicated  only  an  additional  pressure  of  TOp.  A reading 
taken  4 minutes  later  showed  it  to  be  only  40^.  A possible  expleination  is 
that  the  water  vapoxir  was  adsorbed  on  the  shock  tube  walls  thereby  lowering 
the  pressure.  It  is  seen  that  a small  amount  of  water  vapour  added  to  the 
argon  test  gas  can  stabilize  the  shock  wave  and  eliminate  all  disturbances. 

At  the  same  time,  the  ionization  relaxation  length  is  drastically  reduced 
very  much  like  the  effects  produced  on  adding  ~ 0.5^  H2  (by  pressure) . This 
is  not  surprising,  as  the  water  vapour  dissociates  very  readily  at  the  existing 
tenperatures  (~  20,000K)  and  in  turn  ionized  hydrogen  is  produced. 

From  the  foregoing  experiments  it  appears  that  sufficient  water  as 
an  irpurity  will  stabilize  an  argon  shock  wave  at  high  shock  Mach  nunfcers,  just 
like  hydrogen.  However,  not  enough  water  vapour  may  give  rise  to  new  instabil- 
ities. The  latter  is  still  under  study  (Ref.  4l) . One  might  conjecture  that 
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the  small  amounts  of  water  vapour  always  present  in  a shock  tube  could  give 
rise  to  the  instabilities  observed  in  pure  monatomic  gas  runs.  Since  the 
size  of  the  shock  tvibe  used  in  the  present  work  is  large,  it  has  a much 
bigger  inner  surface  area  con^ared  to  the  tubes  used  by  Alpher  and  White 
(Ref.  36)  and  Bershader  et  al  (Refs.  5,  22).  With  a much  bigger  inner  siirface 
area,  the  amoimt  of  adsorbed  water  vapour  is  correspondingly  increased.  Con- 
sequently, such  distiirbances  are  observed  only  in  the  larger  tube.  By  purposely 
increasing  the  amount  of  water  vapour  addition  they  are  eliminated.  Neverthe- 
less, the  disturbances  shown  in  Figs.  4l  and  42  are  of  a very  different  type 
from  those  shown  in  Figs.  8 and  31-  Such  behaviour  would  have  to  be  explained 
by  collision  theory.  Since  the  data  are  few,  more  will  be  known  when  the  work 
of  Ref.  4l  is  con5)leted. 

5 .4  Wall  Effects  in  Relaxation  Region 

An  interesting  result  was  found  (Ref.  30)  experimentally  in  the 
relaxation  region  near  the  shock-tube  wall  as  shown  very  clearly  in  Fig.  8 
for  Xi  = 694 3A,  which  is  more  sensitive  to  electron  density.  It  is  seen 
that  the  electron  front  E moves  in  towards  the  translational  shock  front  S 
near  the  wall  surface.  Similar  wall  effects  can  also  be  found  when  hydrogen 
(Fig.  31) j helium  and  oxygen  (Fig.  13)  imptirities  were  added  to  the  argon 
test  gas  (Ref.  ll).  The  same  effect  is  also  noticed  in  the"  krypton  runs 
(Figs.  9 and  32).  In  the  cases  with  hydrogen  as  an  inpurity  (Figs.  31  and 
32) , the  relaxation  region  is  shortened  so  that  only  a few  fringes  supply 
data  in  the  wall- relaxation  region,  and  therefore  the  change  is  not  so 
apparent . 


The  reasons  for  this  premature  ionization  close  to  the  wall  are 
not  clear.  One  considered  possibility  was  that  a gas-surface  interaction 
occurred  between  the  argon  plasma  and  the  chromium-plated  steel  shock  tube 
wall.  Two  experiments  were  carried  out  by  Brimelow  (Ref.  30)  to  try  to 
eliminate  this  possibility.  In  one,  a covering  (O.OO8  cm)  of  ordinary 
Sellopac  (cellophane)  tape  was  stuck  to  the  wall  and  in  the  other  a thin 
tujTgsten  foil  (0.003  cm)  was  fixed  to  the  wall.  It  was  hoped  that  the  tape, 
being  a good  electrical  insulator,  or  the  tungsten,  having  a different  work 
function  from  .chromium  or  steel,  might  change  the  electron  distribution,  if 
indeed  a gas-surface  interaction  occurred,  however,  no  change  was  observed. 
The  rough  layer  of  sodium  chloride  crystals  on  the  tube  surface  seems  to  bring 
the  two  fronts  even  closer  (Fig.  30).  It  should  be  pointed  out,  however,  that 
only  changes  were  made  in  the  wall  surface  near  the  test  section.  It  is 
dovibtful  if  changes  in  the  wall  surface  over  longer  lengths  of  channel  would 
be  more  effective.  If  adsorbed  water  vapoior  on  the  shock-tube  walls  cause 
this  phenomenon.  Just  as  it  reduces  the  relaxation  length  in  the  free  stream, 
then  it  will  always  be  present  no  matter  what  the  surface  may  be.  One 
difficulty  arises  from  the  fact  that  near  the  wall  the  viscous  boundary  layer 
is  cold  and  the  degree  of  ionization  is  approaching  zero.  However,  at  the 
leading  edge  right  at  the  translational  shock  front,  although  the  boundary 
thickness  is  approaching  zero  and  the  translational  tenperatirre  is  even  higher 
than  at  the  electron  cascade  front,  the  electron  density  is  approaching  zero. 
Consequently,  it  is  difficult  to  give  a substantive  reason  for  the  approach 
of  the  electron  cascade  front  towards  the  translational  shock  front  near  the 
wall. 
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5.5  Sovirces  of  Errors 


In  Ref.  11  there  is  a detailed  discussion  of  the  sources  of  error 
Involved  in  the  present  type  of  experimental  work.  They  can  he  generalized 
as  (l)  errors  associated  in  using  the  interferometric  method  and  (2)  errors 
involved  in  the  measurement  of  the  initial  conditions. 

In  plasma  flows  with  large  density  gradients,  such  as  the  boundary 
layer,  the  light  rays  are  refracted  and  travel  over  longer  optical  paths 
where  the  plasma  refractive  index  is  constantly  changing.  This  introduces 
errors  that  must  be  accounted  for.  In  the  shock  structure  and  in  the  induced 
quasi-uniform  flow  such  gradients  are  considered  as  small.  The  test  section 
side  window  boundary  layer  effects  are  also  neglected  in  the  analysis.  Cal- 
culation done  by  Bunting  and  DeVoto  (Ref.  h2)  on  refractive  errors  involved 
in  interferometric  st\idies  showed  that  the  above  mentioned  errors  were  rela- 
tively smaller  than  the  fringe  shift  reading  error  (typically  ± 0.08  in  the 
present  study) . It  is  therefore  reasonable  to  neglect  such  errors  in  the 
data  reduction.  The  fact  that  in  the  case  of  a plasma,  electrons  and  heavy 
particles  have  opposite  effects  on  light  refraction,  further  reduces  possible 
errors  from  this  source. 

Since  the  locations  of  fringes  are  determined  visually,  it  is  con- 
sidered to  be  a major  error  source  in  the  data  reduction.  The  more  the 
iirtcrferograms  are  enlarged,  the  more  accurate  are  the  fringe  readings. 

However,  since  the  entire  interferogram  is  evaluated,  it  is  inpractical  to 
blow  it  up  to  an  excessive  size.  Normally  it  was  enlarged  to  about  50  cm  x 
50  cm.  A fringe  error  of  ± 0.08,  which  is  considered  to  be  generous,  was 
then  used  to  determine  the  error  associated  ^th  the  flow  properties.  The 
error,  in  the  argon  cases,  was  about  3 x 10“°  gm/cm3  or  about  3.5^  of  the 
deduced  maximimi  density  and  5 x or  jfo  of  the  maximiim  electron  number 

density  when  the  initial  test  gas  pressiire  was  5 torr.  Since  the  above  errors 
depend  on  the  initial  conditions,  they  are  uniform  over  the  entire  flow.  The 
error  in  the  degree  of  ionization,  however,  depends  on  both  local  density  and 
electron  number  density,  so  it  varies  at  different  locations  and  is  usually 
about  0.4%  or  less.  The  errors  due  to  fringe  readings  in  the  krypton  cases 
are  foxuid  to  have  about  the  same  magnitude  as  those  for  argon. 

The  maximum  reading  error  in  the  initial  terperature  was  ± 0.5C, 
and  has  practically  no  effect  on  the  calculated  post-shock  properties. 

The  measurements  of  the  initial  pressure  were  fotind  to  have  mainly 
visual  reading  errors.  A maximvim  random  error  of  about  1%  in  the  pressure 
reading  was  possible.  The  higher  the  initial  presstire,  the  lower  the  error. 
Igra  (Ref.  38)  showed  that  when  the  maximvim  possible  error  in  initial  pressure 
was  considered,  all  calculated  post- shock  properties  were  within  5%  of  the 
nominal  values. 

The  shock  arrival  was  monitored  with  Atlantic-Re search  LD-25  pressure 
transducers  which  triggered  a microsecond  timer.  Considering  both  the  rise 
time  of  the  transducers  and  the  timer  resolution,  the  error  in  the  measured 
time  interval  is  ± 1 jisec.  Depending  on  the  shock  velocity,  the  Mach  number 
so  determined  is  within  ±0.3.  Section  5*1  discussed  the  effect  of  Mach 
nvmber  on  the  deduced  atom-atcm  excitation  collision  cross  section.  The 
effects  of  the  Mach  number  error  on  the  calculated  post- shock  properties  can 
be  seen  in  Figs.  15  and  I6,  and  it  is  significant. 
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Another  type  of  error  in  the  data  reduction  arises  from  neglecting 
the  refractive- index  effects  of  the  added  hydrogen  inpiarities.  However,  in 
view  of  the  fair  to  good  agreement  between  the  analysis  and  the  measured  data 
in  the  cases  involving  the  small  amounts  of  hydrogen  inpurity  (Figs.  33-^0), 
this  error  is  believed  to  have  very  little  effect. 


6.  CONCl.USIONS 


The  use  of  two-wavelength  interferometry  has  proven  to  be  a very 
effective  method  for  investigating  ionizing  shock  wave  structure  in  monatomic 
gases.  The  measurements  of  total  density  and  electron  number  density  in  the 
relaxation  region  provided  the  data  for  evaluating  the  atom-atom  excitation 
collision  cross-sections  for  argon  and  krypton  as  1.0  x 10”19  an^/eV  and 
1.2  X 10-19  cm^/eV,  respectively. 

The  measurements  in  the  radiative  cooling  region  also  provided 
the  information  required  to  choose  an  appropriate  theoretical  radiation 
model.  It  is  shown  that  radiative  cooling  losses  are  important  for  high 
Mach  number  nans  at  higher  pressoires.  The  particle  velocity,  the  electron 
temperature  and  the  atom  temperature  decrease  monatonically  throughooit  the 
radiative  cooling  region,  whereas  the  density  and  pressure  increase.  The 
overall  agreement  between  theory  and  experiment  is  good. 

Oxygen,  helium  and  sodium  chloride  impurities  when  added  to  argon 
do  not  seem  to  have  any  effect  on  shock-wave  structure  and  stability. 

Owing  to  the  low  mass  of  the  hydrogen  atom  and  the  large  excitation 
cross-section  constants  between  the  hydrogen  atom-atom  and  the  hy'drogen 
atom-electron  collisions,  even  a small  addition  (0.i+^)  of  hydrogen  to  the  test 
gas  drastically  reduces  the  relaxation  length.  The  addition  of  hydrogen  also 
removes  the  sinusoidal-type  instabilities  from  the  shock  front  as  well  as 
those  from  the  ionization  front  and  the  subsequent  freestream.  How  hydrogen 
removes  these  instabilities  has  as  yet  not  been  explained,  despite  the  accxorate 
analytical  predictions  of  density,  electron  concentration,  and  degree  of  ioniza- 
tion profiles.  It  should  be  noted  that  the  present  experiments  are  not  in  a 
Mach  nvmiber  range  where  real  gas  effects  cause  the  second  derivative  of  pressure 
with  respect  to  volume,  at  constant  entropy  (S^P/^v^)g  to  became  negative,  in 
violation  of  the  second  law  of  themnodynamics,  which  may  give  rise  to  shock- 
wave  instabilities  (Refs.  ik-Vj) . Consequently,  it  must  be  concluded  that 
the  instabilities  encoointered  in  the  present  and  earlier  research  are  associated 
with  a significant  degree  of  ionization.  This  phenomenon  is  still  being  inves- 
tigated. 


Water  vapour  as  an  added  impoority  has  some  very  unpredictable  effects 
on  the  structiare  and  stability  of  shock  waves.  While  it  can  stabilize  dis- 
turbances and  reduce  the  relaxation  time  like  hydrogen,  it  can  also  produce 
new  fonns  of  distoarbances. 

Another  unsolved  problem  is  the  progression  of  the  cascading  electron 
(ionization)  front  towards  the  translational  shock  front  as  the  shock  tube  wall 
is  approached.  The  present  interferometric  profiles  can  be  considered  as  very 
accurate  experimental  data.  As  such,  it  should  stimoolate  additional  theoretical 
research  on  shock- wave  stability  and  premature  ionization  near  shock  tube  walls. 
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FIGURE  1 DEFLAGRATING  COMBUSTION  PRESSURE  TRACE 


^mixture  * ^*9  x in'*  Torr  (15%  H2  ♦ 7.5%  O2  ♦ 72.5%  He  ♦ 5%  H2) 

'■burst  ■ ^ 

Oscilloscope  Setting:  Vert;  1 volt/di v. 

Horiz:  5 msec/div. 

Charge  Amplifier  Setting  : Transducer  Sensitivity  ; 1.09 

Range  : 1,000  psi/volt 


FIGURE  2 SCHEMATIC  OF  MACH-ZEHNOER  INTERFEROMETER 
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FIGURE  4 FLOW  PRESSURES  AT  THE  CENTRE  OF  THE  TEST-SECTION 

a)  Ar;  Pi  - 5.14  Torr,  T^  - 293.6  K,  M « 15.9, 

Oscilloscope  Setting:  Vert.;  1 volt/div. , Horiz.:  20  psec/div, , 
Charge  Amplifier  Setting:  Sensitivity:  1.09,  Range:  10  psig/volt. 

b)  Kr:  Pi  - 5.15  Torr,  Tj  - 298.3  K,  M - 15.0 

Oscilloscope  Setting:  Vert.;  1 volt/div.,  Horiz.:  50  psec/div.. 
Charge  Amplifier  Setting:  Sensitivity;  1.09,  Range:  10  psig/volt. 

S • Translational  and  lonizational  Shock  Front 
C • Contact  Front 


FIGURE  5 


POLAROID  SHOT  OF  NaCl  CRYSTAL  COATING 
ON  THE  TEST-SECTION  LOWER  WALL  PRIOR  TO  FIRING 


PfC-SNOCK  REONMN 


FIGUPP.  6 

TYPICAL  FLOW-INTERFEROGRAM  SHOWING  PRE-SHOCK  REGION,  TRANSLATIONAL  SHUCK  FRONT  S, 
IONIZATION  RELAXATION  REGION  Xg,  ELECTRON  CASCADE  FRONT  E AND  POST- SHOCK  REGION. 


FIGURE  7 «)  FLOW-INTERFEROGRAM  TO  BE  DIGITIZED  SHOWING  LOCATIONS  OF  CROSS-HAIRS 


FIGURE  7 b)  DIGITIZED  VERSION  OF  IFfTERFEROGRAM  SHOWING  ORIGIN  OF  (x,y)  COORDINATE  SYSTEM  X mil) 


jjjtt-WJatOGRAM  OF  A SHOCK  WAVE  Mg  = 15-2  TRAVELLIMG  INTO  KRYPTON  (Pi  = 5.07  TOHR,  T^ 
295.it  K)  FROM  RIGHT  TO  LEFT.  SINUSOIDAL  TYPE  SHOCK  FROKE  AND  INSTABILITIES  IN  THE 
ELECTRON  CASCADE  FRONT  AND  DOWNSTREAM  OF  IT  ARE  ALSO  SEEN.  S DENOTES  TRANSLATIONAL 
SHOCK  FRONT  AND  E DENOTES  EIECTRON  CASCADE  FRONT . 


FIG.  10  (a)  VARIATIONS  OF  HEAVY  PARTICLE  TEMPERATURE  Ta*  ELECTRON  TEMPERATURE  Te 
AND  HEAVY  PARTICLE  VELOCITY  u WITH  DISTANCE  x BEHIND  A SHOCK  FRONT 
(FOR  KRYPTON  TEST  GAS,  CASE  7).  Mg  = 15.1,  Pi  = 5.15  TORE,  T^  = 
296.2  K AND  ui  = 3.33  X 105  ct^sEC. 


FIG.  10  (b)  VARIATIONS  OP  HEAVY  PARTICLE  TEMPERATURE  Ta,  ELECTRON  TEMPERATURE  Te 
AND  HEAVY  PARTICLE  VELOCITY  u WITH  DISTANCE  x BEHIND  A SHOCK  FRONT 
(FOR  ARGON  TEST  GAS,  CASE  l) . Mg  = 15.9,  Pi  = 5.l4  TORR,  Ti  = 293.6  K 
AND  ui  = 5.08  X 105  caVSEC. 


FIG.  12  miEKFEROGRAM  OF  A STRONG  SHOCK  WAVE  Mg  = l6.1  TRAVELLING  FROM  RIGHT 
TO  LEFT  INTO  ARGON  TEST  GAS  (Pl  = 5-15  TQRR  AND  Ti  = 295.9  K) . 
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FIG.  l4  lOTEElFEROGRAM  OF  A SHXK  WAVE  Mg  = 15-9  TRAVELLIUG  FROM  RIGHT  TO  LEFT 
INTO  ARGON,  = 5-14  TORR  AND  Ti  = 293.6  K (CASE  l) . 


FIG.  15  VARIATIONS  OF  NONDIMENSIONALIZED  ELECTRON  NUMBER  DENSITY  (ng/ne  e)  ^ 
DENSITY  (p/pe)  WITH  DISTANCE  x FOR  CASE  1.  ’ 

— Ms  = 16.2  AND  = 0.9  X 10-15  cs^/qy 

Ms  = 15.9  and  = 1.0  X 10"  19  can^/eV 

Mg  = 15.6  AND  = 1.1  X 10-19  cm2/eV 

Pi  = 5 -14  TQRR  AND  Ti  = 293.6  K. 


Ms  « 16  2 

Ms  « 15.9 

Ms  » 15,6 


1 


Cose  1 


X(cm) 

FIG.  16  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  (a/Op)  WITH  DISTANCE 
X FOR  CASE  1. 

— Mg  = 16.2  AND  S^_^  = 0.9  X 10‘^9 

Mg  = 15.9  AND  Sj^_^  = 1.0  X 10’^9  om^/eV 

Mg  = 15.6  AND  = 1.1  X 10"^^  an^/eV 

P,  = 5-14  TORR  AND  T^  = 293.6  K. 


PIG.  19  VARIATIONS  OF  NONDIMENSIONALIZED  ELECTRON  NUMBER  DENSITY  (n  /n  ) 
AND  DENSm  (p/p„)  WITH  DISTANCE  X FOR  CASE  3.  e e,iL 

ti 


FIG.  23  VARIATIONS  OF  NONDIMENSIONALIZED  ELECTRON  NUMBER  DENSITY  {%/%  g) 
AND  DENSITY  (p/Pg)  WITH  DISTANCE  x FOR  CASE  5-  ’ 


X (mm) 


FIG.  24  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  {oc/°0  WITH 
DISTANCE  X FOR  CASE  5 . 


FIG.  25  DCTERFKkOGRAM  OF  A SHOCK  WAVE  Mg  = 15.2  TRAVELLING  FROM  RIGHT  TO 
LEFT  INTO  KRYPTON,  Pi  = 5-07  TORR  AND  Ti  = 295.4  K (CASE  6) . 


FK.  26  VARIATIONS  OF  NONDBENSIONALIZED  ELECTRON  NUMBER  DENSITY  (n  /n  „) 
AND  DENSITY  (p/pg)  WITH  DISTANCE  x FOR  CASE  6,  e e,E 


Cose  6 


X(mm) 


Case  6 


FIG.  27  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  (o/oe)  WITH 
DISTANCE  X FCR  CASE  6 


FIG.  28  VARIATIONS  OF  NONDIMEWSIONALIZED  ELECTRON  NUMBER  DENSITY  ( n /n  ) 
AND  DENSITY  (p/pg)  WITH  DISTANCE  x FCK  CASE  7.  e ejE' 


FIG.  29  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  ((Vcto)  WTCH 
DISTANCE  X FOR  CASE  7- 


FIG.  30  IMEERFEKOGRAM  OF  A SHOCK  WAVE  Mg  = 15.7  TRAVELLING  FROM  RIGHT  TO  LEFT 
INTO  NaCl  COMTAMINATED  ARGON,  Pi  = 5-1^  TORR  AND  Ti  = 297.^  K. 


(b)  Ms  = 15.0,  Pi  = 5.21  TORE  (O.38fo  H2)  AM)  Ti  = 295.4  K. 

(c)  Mg  = 15.3,  Pi  = 5.19  tore  (0.6^  H2)  AM)  Ti  = 296.6  K. 

(d)  Mg  = 15.0,  Pi  = 5.13  tore  (1.0^  Hg)  AMD  Ti  = 297-7  K. 


FIG.  33  VARIATIONS  OF  NOHDIMEaJSIONALIZED  ELECTRON  NUICER  DENSITY  (n  /n  ) 
AND  DENSITY  (p/pg)  WITH  DISTANCE  x FOR  CASE  8.  e e,E 


8 12  1.6  2.0  2. 


6 20  2 


FIG.  3^  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  (a/cfc.)  WITH 
DISTANCE  X FOR  CASE  8. 


FIG.  35  VARIATIONS  OF  NONDIMENSIONALIZED  ELEICTRON  HUMBER  DENSITY  (n  /n  ) 
AND  DENSm  (p/pg)  WITH  DISTANCE  x FOR  CASE  9.  e e,E 


Case  9 


X (cm) 


Case  9 


FIG.  36  VARIATION  OF  NONDBEHSIONALIZED  DEGREE  OF  IONIZATION  (o/ofe) 
WITH  DISTANCE  X FOR  CASE  9. 


FIG.  4l  INTERFEROGRAM  OF  A SHOCK  WAVE  TRAVELLING  FROM  RIGHT  TO  LEFT  INTO  ARGON 
(Pl  =5.11  TORR  AND  Tx  = 296.2  K) . TEST  GAS  WAS  CONTAMINATED  BY  A 
"SMALL"  BUT  UNKNOWN  AM)UNT  OF  WATER  VAPOUR.  MACH  NUMBER  Mg  RECORDED 
WAS  18.5,  BUT  IS  BELIEVED  TO  BE  ERRONEOUS  PERHAPS  DUE  TO  OSCILLATING 
SHOCK  FRONT . 


FIG.  42  IMTERFEROGRAM  OF  A SHOCK  WAVE  TRAVELLING  FROM  RIGHT  TO  LEFT  INTO  PURE  ARGON 


FIG.  43  lOTERFEROGRAM  OF  A SHOCK  WAVE  Mg  = l4.9  TRAVELLING  FROM  RIGHT  TO  LEFT  INTO 
ARGON  WITH  WATER  VAPOUR  IMPURITY,  Pi  = 5.12  TORR  AJ©  Ti  = 296. 0 K.  HERE 
0.25  ml  OF  DISTILLED  WATER  WAS  INTRODUCED  INTO  THE  TUBE,  AND  PRESSURE  GAUGE 
REGISTERED  30^l  OR  0.1%  WATER  VAPOUR  BY  PRESSURE. 


TABLE  1,  GAS  PRESSURE  VARIATION  WITH  TIME 


TABLE  3 INITIAL  CONDITIONS  FOR  CASES  STUDIED 


1 

Interferograms 

from 

Present 

work 

Present 

work 

Brimelow, 

Ref. 

30 

Whitten, 

Ref. 

29 

Bristow, 

Ref. 

11 

Present 

work 

Present 

work 

Brimelow, 

Ref. 

30 

Present 

work 

Present 

work 

Present 

work 

Ms 

Pj  (Torr) 

i 

Tj  (K) 

Cases 

15.9 

5.14 

293.6 

Ar 

16.1 

5.15 

295.9 

Ar 

16.5 

5.12 

296.6 

Ar 

13.0 

5.01 

296.6 

Ar 

17.2 

2.60 

293.3 

Ar 

15.2 

5.07 

295.4 

Kr 

15.1 

5.15 

296.2 

Kr 

16.2 

5.17 

297.0 

Ar  + 0.4%  H2 

17.7 

3.12 

298.2 

Ar  ♦ 0.4% 

15.0 

5.21 

295.4 

Kr  0.38%  H2 

15.1 

5.16 

297.1 

Kr  ♦ 0.38%  H2 

TABLE  4 

THERMAL  EQUILIBRIUM  VALUES  OF  THE  DENSITY  Pg, 

ELECTRON  NUMBER  DENSITY  AND  raCPEE  OF  IONIZATION  Og  WITHOUT  RADIATION  LOSSES 


APPEMDIX  A - ADDITIONAL  RESULTS 


The  results  presented  in  this  Appendix  were  obtained  from  the  same 
interferograms  (but  at  different  locations)  used  in  Case  1 to  Case  11.  The 
initial  conditions  and  thermal  equilibrium  values  can  be  found  in  Tables  3 
and  4.  These  results  have  been  discussed  in  Section  5- 


A-1 


FIG.  A1  VARIATIONS  OF  NONDIMENSIONALIZED  ELEX3TRON  NUMBER  DENSITY  e)  ^ 

DENSITY  (p/pe)  VOTH  DISTANCE  x FOR  CASE  1.  ’ 

— Ms  = 16.2  AND  sJp.Ar  = 0-9  x 10-19  cn^/ev 

Mg  = 15.9  and  S^.Ar  = 1-0  X IQ-l^  a^/eV 

Mg  = 15.6  AND  = 1.1  X 10-19  cB^/eV 

Pi  = 5-1^  TQRR  and  Ti  = 293.6  K. 


X (mm) 

FIG.  A2  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  (ot/oftr)  WITH  DISTANCE 
X FOR  CASE  1. 

— Mg  = 16.2  AND  = 0.9  X 10"^9 

Mg  = 15.9  and  = 1.0  X 10'^^  cm^/eV 

-•-Mg  » 15.6  AND  = 1.1  X 10"^9  cm^/eV 

Pj^  = 5.1^  TORR  AND  - 293.6  K. 


FIG.  A3  VARIATIONS  OF  NONDIME3rSIONALIZE3)  ELECTRON  NU1«ER  DENSITY  (ne/°e  e)  ^ 
DENSITY  (p/pe)  WITH  DISTANCE  X FOR  CASE  1.  ’ 

— Ms  = 16.2  AND  sJr-Ar  = 0*9  x 10*19  cm^/eV 

Mg  = 15.9  and  S^_Ar  = 1-0  X 10*19  cn^/eV 

Mg  = 15.6  AND  = 1.1  X 10-19  cm^/eV 

Pi  = 5-1^  TCRR  and  Ti  = 293.6  K. 


FIG.  Ah  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  (o/Cte)  WITH  DISTANCE 
X FOR  CASE  1. 

— - 16.2  AND  = 0.9  X 10"^9  cm^/eV 

Mg  - 15.9  AND  = 1.0.  JC  10"^^  cm^/eV 

M.  - 15.6  AND  St.  - 1.1  X 10"^9  ctt^/eV 

- 5.1**  TQRR  and  Tj^  » 293 .6  k. 


FIG.  A7  VARIATIONS  OF  NONDIMEMSIONALIZED  ELECTRON  NUMBER  DENSITY 

AND  DENSITY  (p/Pg)  >nTH  DISTANCE  x FOR  CASE  3*  ’ 


oScf  ^ lOmZATION  (q/<^) 


fig.  A13  VARIATIONS  OP  NOHDIMENSIONALIZED  ELECTRON  NUMBER  DENSITY  (n  /n  „) 
AND  DENSITY  (p/Pg)  WITH  DISTANCE  x FOR  CASE  6.  e e,E 


Case  6 


X (mm) 


30  40 

X ( mm ) 


PIG.  Ali)  VARIATION  OF  NONDIMEHSIONALIZED  DEGREE  OF  IONIZATION  (o/ote)  WITH 
DISTANCE  X FOR  CASE  6 


FIG.  A27  VARIAII0M5  OF  NONDBffiNSIONALIZED  EIECTRON  NUMBER  DENSITY  (n  /n  ) 
AND  DENSITY  (p/pg)  WITH  DISTANCE  x FOR  CASE  10.  e- 


FIG.  A28  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  (o/cL.)  WIOI 
DISTANCE  X FCR  CASE  10.  ^ 


FIG.  A29  VARIATIONS  OF  NONDIMEajSIONALIZED  ELEX3TRON  NUMBER  DENSITY  (n  /n  _) 
AND  DENSITY  (p/Pg)  WITH  DISTANCE  x FOR  CASE  10.  ^ 


X (mm) 

FIG.  A30  VARIATION  OF  NONDIMENSIONALIZED  DECIREE  CF  IONIZATION  (o/gu)  WITH 
DISTANCE  X FOR  CASE  10.  ‘ 


A33 


VARIATIONS  OF  NONDIMENSIONALIZED  ELECTRON  NUMBER  DENSITY  (n  /n  ) 
AND  DENSm  (p/Pg)  VaiH  DISTANCE  x FCR  CASE  11.  S'  e,B 


Case  II 


FIG.  A35  VARIATIONS  OF  NONDIMESiSIONALIZED  ELECTRON  NUMBER  DENSIK  (n  /n  _) 
AND  DENSITY  (p/Pg)  WITH  DISTANCE  x FOR  CASE  11.  e e,JSi 


X (mrr) 


FIG.  A36  VARIATION  OF  NONDIMENSIONALIZED  DEGREE  OF  IONIZATION  (o/cl)  WITH 
DISTANCE  X FOR  CASE  11. 
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entire  flow  and  at  the  same  time  drastically  chorions  the  relaxation  length  to  about  1/3  of 
the  original  value  in  orgon  and  2/3  in  krypton.  The  addition  of  about  1%  of  water  vapour 
(by  pressure)  to  the  argon  test  gas  also  stabilizes  the  flow  and  shortens  the  relaxation 
region.  A thin  coating  of  dissolved  sodium  chloride  in  water  on  the  shock  tube  wall, 
which  on  evacuation  leaves  a rough  surface  of  crystals,  however,  does  not  seem  to  have 
any  effect  on  the  flow.  Although  the  small  impurities  of  hydrogen  and  water  vapour  have 
a significant  effect  in  stabilizing  the  shock,  the  overall  values  of  density  and  electron 
number  density  in  the  quasi-equilibrium  region  remain  nearly  the  same.  The  present 
experimental  results  provided  the  data  for  a theoretical  study  which  took  into  account 
atom-atom  collisions,  electron-atom  collisions  and  radiation-energy  losses  to  deduce  a 
precise  value  for  the  atom-atom  collision  excitation  cross-section  constants  for  argon  and 
krj^ton  having  values  of  S*Ar-Ar  =1.0  x 10-19  cm^/eV  and  S*](i.-Kr  =1*2  x 10“19 
cmVeV,  respectively. 


